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"Measurement of Overload Effect in Photomultipliers as 
Applicable in Laser Radar Techniques" 

D. P. Jones and G. S. Kent 
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University of the West Indies 

"Laser Doppler Measurement of Atmospheric Wind 

Velocity" R. L. Schwiesow, N. L. Abshire, 

and V. E. Derr 
Wave Propagation Laboratory 
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National Oceanic and Atmospheric Administration 

Boulder, Colorado 
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Chairman, B. M. Herman 
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Matrix Elements of Polydisperse Systems of Irregular 
Particles With the Matrix Elements Calculated for 
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N ASA A/Val lops Station 
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P. J. Hof man and W. A. Munn 
Smithsonian Astrophysical Observatory 

"Determination of Atmospheric Aerosol Characteristics 
From the Polarization of Scattered Radiation" 
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Old Dominion University 
and 
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"Inversion of Actual Angular Scattering Data to Obtain 
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Gerald W. Grams 
National Center for Atmospheric Research 
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"Technique for Obtaining Vertical Profiles of Backscatter- 
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J. D. Spinhirne, B. M. Herman , 
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"An Approximate Equation for the Multiply Scattered 
Contribution to a Lidar Return" Edwin W. Eloranta 
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University of Wisconsin 
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M. L. Wright and L S. Gasiorek 
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SESSION III - Stratospheric Scattering Measurements 

Chairman, B. G. Schuster 

3:30- 3:45 Invited Paper, "The Global Stratospheric Aerosol Burden" 28 

J. M. Rosen , D. J . Hofmann, 
and T. J \ Pepin 
University of Wyoming 

3:50- 4:00 "Lidar Measurements of Stratospheric Particulates" 29 

R. D. Hake, Jr., W. Viezee, 
and P. B. Russell 
Stanford Research Institute 
Menlo Park, California 

4:05- 4:15 "Comparison of Lidar Aerosol Measurements and Direct 31 

Particulate Sampling" Frederick G. Fernald, 

Burton G. Schuster, and Charles L. Frush 
National Center for Atmospheric Research 
Boulder, Colorado 

4:20— 4:30 "Global Reconnaissance of Stratospheric Aerosol by 34 

Airborne Lidar" Burton G. Schuster, 

Frederick G. Fernald , and Charles L. Frush 
National Center for Atmospheric Research 
Boulder, Colorado 

4:35— 4:45 "Comparison of LIDAR and In-Situ Measurements of 35 

Stratospheric Aerosols" S. H. Melfi, 

G. B. Northam, and M. P. McCormick 
NASA Langley Research Center 

and 

J. M. Rosen, T. J. Pepin, and D. J. Hofmann 
University of Wyoming 

4:50— 5:00 "Stratospheric Aerosol Lidar Measurements Over 36 

Jerusalem" A . Cohen and M. Graber 

Department of Atmospheric Sciences 
The Hebrew University of Jerusalem, Israel 

5:05- 5:15 "Laser Radar Observations on the Mesosphere and Lower 37 

Thermosphere Over Kingston, Jamaica" 

G. S. Kent and W. Keeniiside 
Department of Physics 
University of the West Indies 
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TIME SESSION IV — Raman, Resonant, and Fluorescent 

Scattering — Laboratory Chairman, C. M. Penney 

8:30— 8:45 Invited Paper, "Measurement of Resonance Scattering Cross 
Section of Sodium D Lines and Laser Radar Detection 
of Sodium Layer in the Upper Atmosphere by a 
Tunable Dye Laser" M. Jyumonji, T. Kobayasi, and 

H. Inaba 

Research Institute of Electrical Communication 
Tohoku University, Sendai, Japan 

and 

T. A ruga and H. Kamiyama 
Institute of Geophysics, Faculty of Science 
Tohoku University, Sendai, Japan 

8:50— 9:00 "Fluorescence of Ambient Atmospheric Aerosols and Its 
Implications for Remote Detection of Gaseous 
Atmospheric Pollutants" Jerry Gefbwachs 

and Milton Birnbaum 
The Aerospace Corporation 
Electronics Research Laboratory 
Los Angeles, California 

9:05— 9:15 "SO 2 Spectroscopy With a Tunable UV Laser" 

W. W. Morey , C. M. Penney , and M, Lapp 
Optical Physics Branch 
General Electric Co., Corporate Research and 
Development, Schenectady, New York 

9:20- 9:30 "Inelastic Light Scattering Processes" Daniel G. Fouche 

M.l.T. Lincoln Laboratory 
Lexington, Mass. 

and 

Richard K. Chang 
Department of Engineering and Applied Science 

Yale University 
New Haven, Connecticut 

9:35— 9:45 "Atmospheric Temperature From Raman Scattering" 

M. Lapp and C. M. Penney 
Optical Physics Branch 
General Electric Co., Corporate Research and 
Development, Schenectady, New York 

9:50—10:00 "Remote Measurement of Atmospheric Temperatures by 

Raman Lidar" Jack A. Safzman 

and Thom A. Coney 
NASA Lewis Research Center 
Cleveland, Ohio 

Coffee 10:05—10:30 a.m. 

BUSINESS SESSION T0:30 a.m. 

Chairman, GLAP Executive Committee J. A. Cooney 

L unc h 1 1 :30 a.m.-l :00 p.m. 
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1 :35— 1 :45 "An Error Analysis of the Differential Absorption 54 

Technique" Richard M. Schotland 

New York University 

1:50— 2:00 "Laser Radar Study Using Resonance Absorption for 57 

Remote Detection of Air Pollutants" 

Takashi Igarashi 
Radio Research Laboratories 
Ministry of Posts and Telecommunications 
Koganei, Tokyo, Japan 

2:05- 2:15 "Determination of Spatial Distribution of Air Pollution by 60 

Dye Laser Measurement of Differential Absorption of 
Elastic Backscatter" S. A. Ahmed 

and J. S. Gergeley 
Electrical Engineering Dept. 

The City College of New York 
New York, New York 

2:20- 2:30 "Remote Sensing of Methane Using an Erbium/Yag Laser: 61 

A Feasibility Study" Kenneth O. White 

and Gerald T. Wade 
Atmospheric Sciences Laboratory 
U.S. Army Electronics Command 
White Sands Missile Range, New Mexico 

and 

Stuart A. Schleusener 
Electrical Engineering Department 
New Mexico State University 
Las Cruces, New Mexico 

2:35— 2:45 "On Atmospheric Temperature Measurement by Lidar 63 

Differential Absorption" James R. Mason 

Atmospheric Sciences Laboratory 
U.S. Army Electronics Command 
White Sands Missile Range, New Mexico 

Coffee 2:50—3:05 p.m. 

3:05— 3:15 "Raman Scattering Frpm Atmospheric Nitrogen in the 65 

Stratosphere" M. J. Garvey and G. S. Kent 

Department of Physics 
University of the West Indies 
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3:20— 3:30 "Laser Radar Measurements of Atmospheric Potassium" 

F. Felix , W. Keertliside, and G. S. Kent 
Department of Physics 
University of the West Indies 
and 

M. C. W . Sandford 
Radio and Space Research Station 
Ditton Park, Slough, England 

3:35— 3:45 "Feasibility of Remote Measurement of SO 2 and NO in 

Smoke Stack Plumes" E. K. Proctor , 

M. L Wright , and Graham Black 
Stanford Research Institute 
Menlo Park, California 


3:50— 4:00 "Detection of SO 2 and NO 2 in Stack Plume by Raman 
Scattering and Fluorescence" 

S. Nakahara, K. Ito, and S. I to 
Mitsubishi Electric Corp., Kamakura Works 
Kamakura, Japan 

4:05— 4:15 "LIDAR Observations of Raman Scattering From SO 2 

in a Power Plant Stack Plume" M. L Brumfield , 

S. H. Melfi, and R. W. Storey , Jr. 

4:20- 4:30 "Raman LIDAR Measurements of Aircraft Turbine Engine 
Exhaust Emissions" D . A. Leonard 


4:35- 4:45 "A Method of Monitoring Non-Resonant Raman Lidar 

Returns During Daylight Hours" John Cooney 

Drexel University 
Physics and Atmospheric Science Department 

Social Hour 6:15 p.m. 

Banquet 7:15 p.m. 

Speaker J. D. Lawrence 

Head, Environmental Quality 
Measurement Systems Office 
LRC 
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SESSION VI — Tropospheric Scattering Measurements 

Chairman, R. T. H. Col Us 
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Stanford Research Institute 
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8:50— 9:00 "Lidar Scattering in the Troposphere" A. /. Carswell, 

J. D. Houston , W. R. McNeil , and S. R. Pa! 
Centre for Research in Experimental Space 
Science and Department of Physics 
York University 
Toronto, Canada 


66 


67 


68 


70 


71 

72 


73 


74 


viii 



TIME 


SESSION VI - Concluded 


PAGE 


9:05- 9:15 "Lidar Measurements of the Thermic Structure in the 76 

Troposphere” Ch. Werner 

German Establishment for Aeronautical and 
Space Research (DFVLR) 
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Brookhaven National Laboratory 
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10:00-10:10 "Estimating Vertical Diffusion Coefficients by Lidar" 79 

Walter M. Culkowski and Searle D. Swisher 
Atmospheric Turbulence and Diffusion Lab. 

National Oceanic and Atmospheric Admin. 

Oak Ridge, Tennessee 

10:15-10:25 "Atmospheric Aerosol and Thermal Structure in the 80 

Boundary Layer Over the Los Angeles Basin" 

Dr, Warren B. Johnson 
Atmospheric Sciences Laboratory 
Stanford Research Institute 
Menlo Park, California 

10:30—10:40 "Laser Sounding of Industrial Mist” V. E. Zuev , 81 

B. V. Kaul, and L V. Samokhvalov 
Institute of Atmospheric Optics of the 
Siberian Branch of the U.S.S.R. Academy of Sciences 

Tomsk, U.S.S.R. 

and 

L. S. Ivlev and K. Ya. Kondratiev 
Leningrad State University 
Leningrad, U.S.S.R. 

10:45-10:55 "Tropospheric Transmissivity Measurements Using the 83 

Raman Nitrogen Lidar Technique” 

M. P. McCormick , G. B. Northam, 
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11:00-11:10 "A Comparison of Atmospheric Structure as Observed With 84 

Lidar and Acoustic Sounder Techniques” 

Edward E. Uthe 
Stanford Research Institute 
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and 
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Argonne National Laboratories 
Argonne, Illinois 

11:15—11:25 "Observations of Invisible Plume by SHGed YAG Laser 88 

Radar” H. Tanizaki, F. !to f T. Yakuo, 
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1-10, Nisshin Cho, Fuchu City, 183 Japan 
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T. Ikeno 

Japan Weather Association 
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Invited Paper, "Detection and Estimation of Atmospheric 
Turbidity From the Polarization of Scattered Radiation" 
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Atmospheric Optics and Radiation Laboratory 
Department of Meteorology, University of California 

Los Angeles, California 

"On the Inversion of Light Scattering Data Into 
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Department of Atmospheric Sciences 
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Boulder, Colorado 
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Tomsk, U.S.S.R. 
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LIDAR SYSTEMS - WHERE WE ARE 

hy 

M. P. McCormick 

NASA - Langley Research Center 
ABSTRACT 

A general review of the past ten years of LIDAR system development will 
he presented. Equipment problems and many of their solutions are discussed. 
These include problems with photomultipliers , lasers , interference filters , 
data acquisition, etc. Calibration and alignment techniques are also pre- 
sented. Instrumentation for multiwavelength , Raman, doppler, monostatic and 
bistatic techniques being used by various researchers are described. New 
developments and a gaze into the future are included. 
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MOBILE REMOTE SENSING LABORATORY 


by 

J. J. Grossman, M. Muramoto and J. Kacin 
McDonnell Douglas Astronautics Company-West 

Huntington Beach, California 92647 
ABSTRACT 

The McDonnell Douglas Astronautics Company (MDAC), has "built a mobile re- 
mote sensing laboratory for studying atmospheric optical phenomena and identi- 
fying various molecular gases present in the atmosphere. 

Major components of this laboratory are (A) the laser system, (b) receiv- 
ing telescope, (C) spectrometer and detectors, (D) and the data displaying and 
recording electronics. 

The transmitter /receiving telescope (T/R) assembly is mounted on an air- 
stabilized isolation table which is housed in a temperature-controlled 22-ft 

trailer truck. This system can be positioned in azimuth and elevation to an 

-3 

accuracy of less than 2 x 10 rad. with the aid of voltage readouts. The 
spectrometer and the laser system are light coupled to a transit-mounted T/R 
assembly through a mirror-prism system that permits use of the laser and de- 
tector systems at fixed position inside the light-tight van enclosure. 

One laser is a modified Spacerays Model FD 22 High Powered Q-switched 
Ruby laser. The ruby is 6.63" long by 3/8 " diameter and Q-switched by a KD*P 
Pockels cell. The output of 6943$ laser line is about 2 joules with a pulse 
length of about 30 nsec. The system is also equipped with the an ADP crystal 
frequency doubler producing more than 0.15 joules at 3471$. The maximum repe- 
tition rate is 2 pulses per sec. The cavity, lamps, and ruby rod are tempera- 
ture controlled by a Dunham-Bush Water Chiller, and the rod is kept at the 
constant temperature. The output from the laser is collimated by a simple 
Galilean telescope to give the beam divergence of less than 0.75 mrad. 

The back-scattered signal receiving module consists of (l) 12" diam. 
f/l6 Cassegrainian telescope, (a) a multi-mirror beam directing system, and 
(5) 3/4 m modified Jarrell -Ash Double Spectrometer. This spectrometer is 
equipped with polychromator housing which holds several photomultipler (PM) 
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tubes at preselected Raman spectral lines. Also provision has been made so 
that the return signals can be PM detected through dichroic filters and 
various temperature-controlled narrow-band pass filters. 

The display and recording consists of (l) the Tektronix Model R k^k and 
Polaroid camera for recording the output of the photodiode laser energy moni- 
tor 5 (2) a Tektronix model R55 6 Dual-Beam Oscilloscope with Type lAk four- 
channel amplifier and Type 1A1 dual-trace plug-in unit and a Coleman 
Oscillatron KD-5 Camera System for recording the output from the photomulti- 
plier tubes. The output of the photo-diode monitor is the trigger start 
pulse for both oscilloscopes. Additionally * three channels presently are 
monitored by preamplifier/Biomation 6l0 B combinations and the data stored on 
magnetic tape. 



A NEW GaAs LASER RADAR FOR ATMOSPHERIC MEASUREMENTS 


by 

R. T. Brown, Sperry Research Center, Sudbury, Mass. 

and 

A. P. Stoliar, Sperry Gyroscope, Great Neck, New York 

ABSTRACT 

A special GaAs lidar using fiber coupled diode lasers was constructed for 
the purpose of measuring the extinction coefficient distribution within a 
large atmospheric volume at a rate compatible with atmospheric kinematics. The 
technique is based on taking backscatter signature ratios over spatial incre- 
ments after the returns are normalized by pulse integration. Essential aspects 
of the lidar design are beam pulse power, repetition rate, detection system 
dynamic range and decay linearity. It was necessary to preclude the 
possibility of eye hazard under any operating conditions, including directly 
viewing the emitting aperture at close distance with a night-adapted eye. 

The electronic signal processing and control circuits were built to allow * 
versatile operations. 

Extinction coefficient measurements were made in fog and clouds using a 
low-power laboratory version of the lidar, demonstrating feasibility. Data 
are presented showing range squared corrected backscatter profiles converted 
to extinction coefficient profiles, temporal signal fluctuations, and solar 
induced background noise. These results aided in the design of the lidar which 
is described. Functional tests of this lidar and the implications relevant 
to the design of a prototype model are discussed. 

This work was jointly sponsored by Sperry Rand Corporation under its 
Independent Research and Development program; the Air Force Avionics 
Laboratory, Wright Field, Dayton, Ohio; and the Naval Ammunition Depot, Crane, 
Indiana . 


; U 



DESIGN AND TESTING OF AN ERBIUM LASER RANGEFINDER 
FOR U&E AS A CEILOMETER 

"by 

Eugene Y. Moroz 

Air Force Cambridge Research Laboratories 
Bedford ? Massachusetts 
and 

Joseph P. Segre and Norman R. Truscott 
American Optical Corporation 
Framingham, Massachusetts 

ABSTRACT 

The use of pulsed laser rangefinders to determine cloud base heights for 
airfield operations have been investigated in the past. Both ruby and neody- 
mium systems have been shown to be excellent laser ceilometers. However, 
these systems present a potential eye hazard even at a considerable distance 
from the source and, therefore, their routine use at airfields is question- 
able . 

This paper discusses an erbium doped glass laser rangefinder that was 
constructed as an experimental laser ceilometer. The laser operates at 1.5^- 
micrometers (pm) and at this wavelength, the output power density is an order 
of magnitude below the proposed Z136 ANSI standard. On this basis the system 
is considered completely eye safe. 

The laser uses a clad silicate glass rod (AO 2101-2 699 ) 5 ^ mm core and 
6 mm cladding diameter, 150 mm long which is pumped by two h x 6 mm flash- 
lamps. It is Q-switched by means of a sapphire flat, rotating at 2^,000 rpm. 
The laser produces a 35 nanosecond (ns) pulse having a peak power output of 
one megawatt . It can be fired automatically at a rate of one pulse per min- 
ute . 

The receiver optics consists of a 12 inch diameter, f:2, reflecting 
paraboloid. The detector is a germanium avalanche photodiode (TI XL 7 6). 

The receiver field of view is 0.1+ milliradian, matching the transmitter beam 
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divergence. The spectral bandpass is from 1.2 to 1.6 pm. The shor.t wave- 
length cutoff is provided by the receiver window (AO IE glass) and the long 
wavelength cutoff by the photodiode spectral response. Due to the low level 
of sky background in this spectral region, a bright cloud background is of 
the same order as the internal noise of the receiver which has a bO MHz band- 
width. 

The system has been extensively tested on a variety of cloud targets in 
order to evaluate its capability as a laser ceilometer. Cloud heights have 
been measured out to 15,000 feet. A signal to noise ratio of five was ob- 
tained when ranging out to clouds at a slant distance of 33,000 feet. 

Comparative tests were made at AFCRL, Hanscom Field, Bedford, Mass- 
chusetts between this system, a ruby laser ceilometer and a rotating beam 
ceilometer (RBC). The ruby laser ceilometer is manufactured by ASEA of 
Sweden and has a 30 ns , 2 megawatt output. The RBC is a standard AN/GMQ-13 
ceilometer using a U00 foot baseline. Polaroids were made of the three 
oscilloscope displays. Cloud base heights were determined from the traces 
using the criterion of Circular N. All three systems were found to be in 
good agreement. The laser systems showed more detailed information on the 
cloud structure than the RBC. 
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ABSTRACT 


A REAL-TIME LIDAR SIGNAL PROCESSOR 

John Oblanas 5 William Viezee , and R. Collis 
Stanford Research Institute 
Menlo Park, California 

A real-time, wide band analog signal processor has been developed at SRI 
for use in lidar studies of the lower atmosphere. 

Unlike previous designs, this processor does not employ time-varied gain 
of the photomultiplier, thus eliminating a potential source of error in the 
lidar data. The processor generates an analog waveform proportional to the 
function 2 £n R (where R is slant range). This waveform is summed with the 
logarithmic output of the lidar receiver to compensate for the inverse range 
square dependence inherent in lidar returns from volume targets. 

When displayed in A-scope form (amplitude vs range) , the signature of 
the range-corrected logarithmic signal provides a direct indication of the 
degree of homogeneity of the atmosphere and of its extinction coefficient. 

This follows from the fact that the rate of change of received signal, Pr , 
with range, R, is a function of the extinction coefficient C according to the 
equation. 

d(£n PrR 2 ) _ _ 

— -'’O ——————— — — 0 

2 dR 

Thus where the atmosphere is homogeneous (a condition that will be indicated 
by the rectilinearity of the displayed trace) O' can readily be derived from 
the trace's slope. 

Similarly, the output of this processor is well suited for use in in- 
tensity-modulated two dimensional presentations of lidar data obtained by 
scanning (i.e., RHI, PPI, Height /Distance , or Height/Time displays), since 
the range correction greatly reduces the dynamic range requirements and facili- 
tates the interpretation of the displayed data by inspection. For such dis- 
plays , a further reduction in dynamic range can be achieved by increasing the 

, T 



slope of the ramp function 2 In R to compensate , at least in general terms , 
for the effects of attenuation* 
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MONOSTATIC LIDAR DATA RECORDING AND 
REDUCAT I ON SYSTEM 

J , D. Spinhirne, M. Z. Hansen and J. A. Reagen 
The University of Arizona , Tucson 
ABSTRACT 

A hardware system has been constructed to perform signal compression , 
analog to digital conversion and digital recording of monostatic lidar returns 
with 3 l single shot accuracy. Emphasis has been placed on simplicity and ease 
of design and construction. 

The system has been built using comparatively inexpensive, commercial 
units that can be incorporated into lidar data acquisition systems with a mini- 
mal amount of design and construction effort. Two techniques are used to 
compress the photomultiplier signal. An accurate high-speed analog multiplier 
is used to compensate for the l/r decrease of the signal. The signal is also 
compressed with a gain switching amplifer which utilizes low cost, monolithic, 
programable gain operational amplifiers to provide gain switching. The 
amplifier switches and settles in 400 ns. and has a bandwidth of 3*5 mhz on all 
channels. Analog to digital conversion is accomplished with a Biomat ion 6l0 
transient recorder which has a versatile input and recording capability and 
provides a 6-bit, 256 word output. 

Data is punched onto paper tape. In addition to the lidar signal; the 
readout of a digital energy monitor, the slant path angle, instrument settings 
and time of observation are also recorded. Interfacing between the recorder, 
punch and auxiliary data source is accomplished with two cards of control logic. 
During an observation, numerous shots are taken at several slant path angles. 

The data is fed into a computer for rapid reduction of the data and presenta- 
tion of atmospheric optical profiles. 
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A PROPOSED MULTI-CHANNEL PHOTON COUNTER WITH 
REAL-TIME DATA ACQUISITION FOR LASER RADAR 


by 

Gerald T. Wade 
Teddy Barber 

Atmospheric Sciences Laboratory 
US Army Electronics Command 
White Sands Missile Range, New Mexico 88002 

and 

Robert Armstrong 
Physics Department 
New Mexico State University 
Las Cruces, New Mexico 88001 

ABSTRACT 

A high-speed, multi-channel photon counting system has been designed for 
use with a variety of laser radar experiments. It utilizes two solid state 
buffer memories to replace the large array of counters normally associated 
with such a system and a minicomputer to handle data manipulation and display. 
The speed, versatility, compact size, and low material cost of the system make 
it a reasonable alternative to systems presently in use. System features 
include a 250 MHz photon counting rate, 100 nanosecond minimum sample interval, 
up to 512 channel data accumulation with sub-nanosecond switching times between 
channels, and a maximum laser transmitter repetition rate of 10 KHz. This 
system offers several advantages over presently available systems, such as 
allowing higher laser repetition rates without sacrificing spatial resolution 
or the total number of data channels, and enabling real-time modification of 
the system operating parameters as a function of incoming data. This paper 
describes the parameters influencing such a system and its basic design. 
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MEASUREMENT OF OVERLOAD EFFECT IN PHOTOMULTIPLIERS 
AS APPLICABLE IN LASER RADAR TECHNIQUES 


by 

D. P. Jones and G. S. Kent 

Department of Physics, University of the West Indies 

ABSTRACT 

Several types of overload effect, or spurious noise pulses, are 
known to exist in photomultipliers. This paper presents the analysis of 
experiments to determine the magnitude of the enhancement which occurs up 
to 1 ms after the tube has received an intense light pulse. The enhance- 
ment appears as a decaying elevation in the magnitude of the dark current 
and is of particular significance in laser radar measurements of atmospheric 
backscatter. 

Intense light pulses were produced from both a pulsed light emitting 
diode and a ruby laser radar system. Although the intensity from the light 
emitting diode was barely sufficient to produce a noticeable overload effect , 
it served as a guide to the intensity of light the photomultiplier tube could 
withstand before overload occurred. The laser experiments where the over- 
load pulse had a much greater intensity showed that the magnitude of the dark 
current enhancement was proportional to the intensity of the incident radia- 
tion. 

This result is of particular importance in relation to the measurement 
of atmospheric density in the lower atmosphere where, unless special pre- 
cautions are taken, the receiving photomultiplier tube may be subjected to 
an initial light pulse whose intensity is much greater than the signal which 
is measured a short time later. 
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LASER DOPPLER MEASUREMENT OF ATMOSPHERIC WIND VELOCITY 


R. L. Schwiesow, N. L. Abshire and V. E. Derr 
Wave Propagation Laboratory 
Environmental Research Laboratories 
National Oceanic and Atmospheric Administration 
Boulder, Colorado 80302 

ABSTRACT 

Our presentation consists of two parts: (l) a summary review of laser 

Doppler principles and applications, and (2) operational design and preliminary 
laboratory tests of a CO^ laser system for NOAA applications. 

Laser Doppler velocimeters for short range applications to water tunnels 
and seeded- flow wind tunnels are widely used and reviewed"'". The non-relativ- 
istic expression for the Doppler shift, 

Vv D = (n/2Tr) v • (k g - k Q ) , 

where n is the refractive index of the medium, v the scatterer velocity 

and k g and k Q the scattered and incident wave vectors (2'ir/X) , applies to 

both local measurement and remote measurement of flows . Remote velocity 

measurements with a single-ended system at ranges of tens of meters or more 

severely restrict the usable experimental configurations, however. We will 

2 3 

critically review the atmospheric results of infrared heterodyne , two-beam 

L 5 

visible differential , and visible incoherent detection systems , as well as 

6 

some signal to noise calculations . 

For the measurement of wind in severe storm and boundary layer meteorology 
applications, considerations of background interference, atmospheric refractive 
index fluctuations, eye safety, laser stability and efficiency, and other 
practical concerns are properly discussed. Although a heterodyne, CO^ laser 
system is attractive for these wind measurement purposes, an analysis of 

7 

expected return based on the atmospheric aerosol model of McClatchey, et al. , 
is required to substantiate the adequacy of the technique. 

We will present an outline of the NOAA Doppler lidar equipment and exhibit 
velocity data from a diffuse scattering source at 12 meters. Conference 
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participants are invited to an interactive discussion of the application of 
this type of laser radar to problems such as waterspouts, thunderstorms, and 
helical rolls. 
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A COMPARISON OF THE MEASURED SCATTERING MATRIX ELEMENTS OF POLYDISPERSE 
SYSTEMS OF IRREGULAR PARTICLES WITH THE MATRIX ELEMENTS CALCULATED FOR 
SPHERES THAT HAVE THE SAME DISTRIBUTION OF PROJECTED AREAS 

o + + 

A. C. Holland, P. J. Hofman and W. A. Munn 

ABSTRACT 

The elements of the scattering or Mueller Matrix for three polydisperse sys- 
tems of irregular, randomly -oriented particles have "been measured in absolute 
terms as a function of scattering angle for several visible wavelengths. 

The samples consisted of commercially available silicon dioxide particles that 
fit three distinct lognormal size distributions. 

The measured matrix elements were compared with the matrix elements calculated 
for spheres that had the same refractive index and fitted the same distribu- 
tion of projected areas. Correlations and discrepancies between the two 
sets of matrix elements will be discussed. 


0 NASA/Wallops Station 
+ Smithsonian Astrophysical Observatory 
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DETERMINATION OF ATMOSPHERIC AEROSOL CHARACTERISTICS 
FROM THE POLARIZATION OF SCATTERED RADIATION 


F. S. Harris, Jr. M. P. McCormick 

Old Dominion University NASA Langley Research Center 

ABSTRACT 

Aerosols affect the polarization of radiation in scattering, hence 
measured polarization can be used to infer the nature of the particles. 

Size distribution, particle shape, real and absorption parts of the complex 
refractive index affect the scattering. From Lorenz-Mie calculations of 
the 4- Stokes parameters as a function of scattering angle for various 
wavelengths the following polarization parameters were plotted: total 

intensity, intensity of polarization in plane of observation, intensity 
perpendicular to the plane of observation, polarization ratio, polarization 
(using all ^-Stokes parameters), plane of the polarization ellipse and its 
ellipticity. A six-component log-Gaussian size distribution model was 
used to study the effects of the nature of the polarization due to variations 
in the size distribution and complex refractive index. Though a rigorous 
inversion from measurements of scattering to detailed specification of 
aerosol characteristics is not possible, considerable information about the 
nature of the aerosols can be obtained. Only single scattering from 
aerosols was used in this paper. Also, the background due to Rayleigh gas 
scattering, the reduction of effects as a result of multiple scattering 

and polarization effects of possible ground background (airborne platforms) 

« 

were not included. 

It is convenient to use radiation initially plane-polarized at 1+5° to 
the plane of observation. The total* parallel and perpendicular component 
intensities are meaningful* "but fluctuations in the source* number of 
particles in 'the scattering volume* and other variables may reduce their 
experimental usefulness. The polarization parameters of polarization ratio* 
polarization* inclination angle and ellipticity are ratios and hence not 
influenced by these changes in the intensities. 

The effect of size distribution can be seen by noting the relative 
contributions to the scattering parameters of the various log-Gaussian 
components. Large particles, for example, give a strong forward peak, a 
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rainbow and enhanced l80° hackscattering in the intensities. Small particles, 
though also different in the parallel and perpendicular components, produce 
a symmetric angular distribution as the sizes get smaller. The variation of 
parallel and perpendicular components with angle and the variation in phase 
difference between them give marked polarization effects especially in the 
backward scattering region up to but not including 180° backscattering. 

At scattering angles of greater than about 70° there are marked differences 
produced by the modified gamma and Junge-type power-law size distribution in 
position and maxima of the polarization ratio, polarization, inclination and 
ellipticity. Appreciable difference in the exponents of the power law 
distributions are distinguishable. Obviously, there are differences in the 
backscattered intensities between Junge and modified gamma size distribution. 

An increase in the real part of the refractive index: moves the 

rainbow to larger scattering angles; enhances the l80° backscattering; 
decreases the polarization ratio maximum (around 90°); and increases the 
slope near 0° inclination angle. An increase in absorption on the other 
hand, flattens and reduces the intensity curves outside the near forward 
scattering region, including the rainbow, glory, and l80° backscattering 
effects. Also the polarization ratio maximum increases, and the inclination 
angle makes a more gradual transition across the 0° inclination axis as it 
moves to smaller scattering angles. These effects are much more evident 
for larger particles. The ellipticity, especially at angles of more than 
100° scattering, is also sensitive to changes in the particle character- 
istics. 

It is obvious that these analytical calculations are useful for 
determining aerosol characteristics. However, in dealing with "real-world” 
aerosols, a multidiscipline approach must be used. This includes: 
measurement of the physical and optical characteristics of the particles 
by various sampling techniques; measurement of the il-Stokes parameters as 
a function of wavelength and angle; and analytical calculations based on 
"real" aerosol characteristics . 
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INVERSION OF ACTUAL ANGULAR SCATTERING DATA 
TO OBTAIN PARTICLE SIZE DISTRIBUTION 

12 1 
Philip B. Russell ’ and Gerald W. Grams 

1 t 2 

National Center for Atmospheric Research Stanford Research Institute 

Boulder, Colorado 80312 Menlo Park, California 94025 

ABSTRACT 

Measurements of the intensity of laser light scattered by a volume of at- 
mospheric aerosol as a function of scattering angle have been numerically 
inverted to obtain the size distribution of particles in the scattering volume. 
The numerical results are compared to simultaneous measurements of particle 
size distribution made by direct sampling techniques. The data were taken on 
ambient soil dust aerosols in Big Spring, Texas, as described by Grams et al. 
(1973). The inversion technique is essentially the "smoothed-least-squares" 
matrix method described by Twomey (1963, 19^5 ) and Westwater and Strand (1972); 
however, modifications were made to handle the large measurement errors present 
in the scattering data. In addition, techniques of singular value decomposition 
(Golub and Reinsch, 1970) were used to minimize errors caused by ill-condition- 
ing of matrices. 

In this method of solution a family of acceptable solutions , rather than 
a unique result;! is obtained; ifee mean of all such solutions, together with 
standard deviations at given radius intervals, is then taken to be the best 
estimate of the particle size distribution. An example of a result obtained 
in this manner, in a case where the dust particles appeared to be scattering as 
homogeneous spheres, is shown in Figs. 1 and 2. A synopsis of the solution 
technique, examples of other cases inverted, and an indication of problems 
presented by particle asphericity, unknown index of refraction, etc., will be 
presented. 

Golub, G. , and C. Reinsch (1970) , "Singular Value Decomposition and Linear 
Least Squares Solutions," Num . Mat . l4 , 403-420. 

Grams, G. W. , I. H. Blifford, Jr., D. A. Gillette, and P. B. Russell (1973), 
"Complex Index of Refraction of Airborne Soil Particles and Implications 
for Climatic Change," submitted to J.A.S. , April 1970. 
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Twomey, S. (1963), "On the Numerical Solution of Fredholm Integral Equations 
of the First Kind by Inversion of the Linear System Produced by 
Quadrature," J.A.C.M . 10, 97-101. 

Twomey, S. ( 1965)9 "The Application of Numerical Filtering to the Solution 
of Integral Equations Encountered in Indirect Sensing Measurements," 

J. Franklin Inst . 279 < 95-109* 

Westwater, E. R. and 0. N. Strand (1972), "Inversion Techniques," Chap. 16 
of Remote Sensing of the Troposphere , V. E. Derr, ed, for sale by Supt. of 
Documents, U. S. Government Printing Office, Wash., D. C. 20k02. 


t The National Center for Atmospheric Research is sponsored by the National 
Science Foundation. 
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FIGURE 1 COMPARISON OF PARTICLE SIZE DISTRIBUTIONS 
OBTAINED BY DIRECT SAMPLING AND INVERSION 
OF LIGHT SCATTERING MEASUREMENTS 
( | direct sampling data point with 

estimated error; lognormal curve 

fitted to sampling data points; | inver- 
sion solution with standard deviation, 
computed from measured cross sections in 
Figure 2) 
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FIGURE 2 TOTAL (GAS PLUS PARTICLE) BISTATIC RADAR 
CROSS SECTIONS 

( | polar nephelometer measurement with 

geometric standard deviation; - computed 

value using inversion solution shown in 
Figure 1) 





TECHNIQUE FOR OBTAINING VERTICAL PROFILES OF 
BACKSCATTERING AND EXTINCTION CROSS SECTIONS 
USING SLANT PATH LIDAR MEASUREMENTS 

J. D. Spinhirne, B. M. Herman and J. A. Reagen 

ABSTRACT 

A method is presented for solving for vertical profiles of atmospheric 
particulate extinction and backscattering cross-sections utilizing monostatic 

1 

lidar slant path measurements. The method is an extension of work by Fernald . 
It is shown that the number of assumptions necessary for an iterative solution 
of extinction and backscattering cross sections can be reduced if lidar slant 
path measurements are used to solve directly for optical depths. 

The technique is useful only if sufficiently accurate lidar measurements 
are available. With highly accurate measurements it is also possible to solve 
directly for extinction cross sections without an iterative solution of a 
transcendental equation if the proper reduction scheme is used. The required 
accuracy is discussed and results showing the effect of errors are presented. 


F. G. Fernald, B. M. Herman and J. A. Reagen, Journal of Applied Meteorology , 
Vol . 11, No. 3 (1973) . 
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AN APPROXIMATE EQUATION FOR THE 


MULTIPLY SCATTERED CONTRIBUTION TO A LIDAR RETURN 

Edwin ¥. Eloranta 
Department of Meteorology 
University of Wisconsin 
Madison, Wisconsin 5370 6 


ABSTRACT 


An approximate equation is developed which describes the contribution of 
th 

N order scattering to a lidar return. This development assumes a homogene- 
ous scattering medium characterized hy a scattering function sharply peaked in 
the forward direction and relatively insensitive to angle near the hackscatter 
direction. The derivation includes the effects of finite divergence of the 
transmitted laser beam, finite receiver field of view* finite separation be- 
tween the laser and the receiver and nonparallel system alignment. 

The derivation presented uses small angle approximations to reduce the 
time dependent multiple scattering problem to a time independent form which 
is then solved with techniques previously developed for multiple small angle 
nuclear scattering. This derivation results in the following equation for 
power received due to order small angle scattering: 
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Power received due to order scattering 
Transmitter energy 
Speed of light 

Scattering cross section per unit volume 
Range 

Area of receiver 

Backscatter phase function (see Deirmendjian 5 1969) 

Order of scattering 

Mean .‘square divergence angle of transmitting laser 

Mean square wisth of the forward scattering peak in the phase 
function 

Modified Bessel function of zero order 
Half-angle receiver field of view 

Lateral separation between laser and receiving telescope 

Angular misalignment of the laser and telescope axis (this misalign- 
ment is assumed to be in the plane defined by the telescope mirror 
and the axis of the laser) 
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Equation 1 provides a simple means of estimating multiple scattering 
during operational lidar work. Notice that in the special case of a nearly 
coaxial system geometry where: ij>«l, I o (2^»^u)-l for all u<a<$£. In this 

case the integral in equation 1 reduces to a simple exponential and 'slide 
rule' solutions can be obtained for P^. 

Examples of the multiply scattered lidar return predicted by equation 1 
will be presented and compared to preliminary measurements. 
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AN OPTICAL MODEL FOR 
STRATOSPHERIC AEROSOLS 


J. M. Rosen 
D. J. Hofmann 
T. J. Pepin 
University of Wyoming 


ABSTRACT 


Most remote sensing techniques of stratospheric aerosol rely on an optical 
model to relate the measured quantities with physical parameters such as con- 
centration, size, and composition. By comparing data taken simultaneously 
using remote sensing and in-situ techniques, it is possible to begin the devel- 
opment of a useful relation between the various results. The heart of this 
comparison is the size distribution and composition of the particles. This 
report describes an optical model that is consistent with impactor samples, 
filter samples (measurement of total mass), photoelectric particle counters, 
LIDAR results, and solar extinction measurements. 



THE EFFECTS OF NON-SPHERICAL PARTICLES ON ESTIMATIONS 
OF AEROSOL PROPERTIES 

By 

Robert ¥. L. Thomas 

Wolf Research and Development Corporation 
6801 Kenilworth Avenue 
Riverdale , Maryland 20801 

ABSTRACT 

A Monte Carlo radiative transfer model has been constructed to investigate 
the transmission and scattering properties of an inhomogeneous but layered 
medium. The program has been applied to investigate the multiple scattering 
properties of the Elterman aerosol model assuming 

1. sperical particles with log-normally distributed radii. 

2. an aerosol with the same mass distribution as (l) but whose scattering 
matrix elements have been determined experimentally by Holland, et al. 
at a wavelength of 0.5^6]J. 

The results are discussed in terms of the detection of non-spherical parti- 
cles by lidar experiments and the possible influence of these particles on 
systematic errors in aerosol property estimations. 
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ABSTRACT 


A MODULAR ATMOSPHERIC PROPAGATION PROGRAM (MAPP) 

by 

M. L. Wright and L. S. Gasiorek 

A computer model of the atmosphere has been developed that will solve a 
wide variety of lidar and atmospheric propagation problems. The lidar system 
calculations involve the prediction of path loss and lidar signal return for a 
variety of atmospheric conditions and lidar system configurations. The program 
contains a very large data base of absorption and attenuation coefficients for 
all of the naturally occurring gases and several gaseous pollutants such as 
SO^, NO^ , NO 9 N^O , HCL, etc. The spectral range of this data base extends 
from the ultra-violet (1200 £) through the infra-red. Several aerosol models 
are included in the program to account for the effects of this important 
scattering mechanism. 

The program incorporates a maximum of flexibility in accepting input 
data and providing output data in an easily used form. Profiles for the 
standard atmosphere and for various polluted atmosphere and cloud models are 
stored in the machine. These stored models can be used singly or in combination, 
or be supplemented or replaced by additional profiles provided as input data to 
the program. A variety of output options are available involving either numeri- 
cal data or graphical data; the output quantities can be attenuation, power 
returned or power differences, as functions of range or path position. 

Path loss is computed by dividing the distance into increments, calculating 
the absorption for each increment and summing up the loss over the desired 
path. The absorption for each increment is calculated by multiplying the 
material concentration times the absorption coefficient for that material. 

The absorption coefficient for a particular wavelength and material is 
determined by either table look-up or calculating the coefficient from spectral 
line data. The line profile is calculated as a function of temperature and 
pressure. The concentration of the gas as a function of range is either taken 
from a stored model or an input quantity. 
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By appropriate selection of a peak absorption and a peak transmission 
wavelength 5 the program will calculate a range-resolved return for the 
differential-absorption lidar (DIAL) measurement technique. Multiple 
computer runs may be made to determine the optimum pair of wavelengths in 
using the DIAL technique to monitor the constituent of interest in the 
presence of other interfering gases. 
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THE GLOBAL STRATOSPHERIC 
AEROSOL BURDEN 

J. M. Rosen, D. J. Hofmann, T. J. Pepin 
University of Wyoming 
ABSTRACT 

Since mid-1972 the University of Wyoming f s Atmospheric Physics Group has 
heen conducting bimonthly global surveys of the stratospheric aerosol using 
balloon-borne detectors. A total of ten stations are included in the sounding 
network, of which there are seven in the northern hemisphere and three in the 
southern hemisphere. From the results of these measurements new information 
is now becoming available concerning the temporal variations and spatial 
distribution of aerosols. An effort is being made to assess the characteristics 
of the natural stratospheric aerosol background. However, since data is avail- 
able only over a limited time period, it is still somewhat difficult to 
determine whether or not the present aerosol concentration levels are in a 
steady state condition. 
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LIDAR MEASUREMENTS OF STRATOSPHERIC PARTICULATES 


R. D. Hake, Jr., W. Viezee , and P. B. Russell 
Stanford Research Institute, Menlo Park, California 

ABSTRACT 

Since October 1972, Stanford Research Institute has conducted a program 
of periodic nighttime observations of stratospheric aerosols using a ground- 
base lidar. The lidar system consists of a 2 J, Q-switched ruby laser 
fired at a rate of 0.5 Hz, and a l6-inch Cassegrain receiver; with 
simultaneous analog and pulse-counting modes of data acquisition. The 
analog data is usable to heights near 30 km, and the pulse-counting data 
is used above 25 km. Typically 20 minutes of operation is required to 
obtain one percent accuracy in 1/4 -km range bins at 25-km altitude. 

The basic data are in the form of vertical profiles of atmospheric 
backscatter in relative units. Techniques to deduce from these data 
information pertinent to determining the spatial and temporal distribution 
of particulate matter in the stratosphere are illustrated. 

Observations in the late fall of 1972 and spring of 1973 are analyzed 
in terms of vertical profiles of the aerosol backseattering coefficient 
(see Figure l) that show the presence and height of stratospheric aerosol 
scattering layers but not their nature. Observed maximum values of total 
(aerosol plus molecular) backscatter are only 10 to 15$ above the molecular 
backscatter. Multiple-layered structure is evident in the vertical 
aerosol distribution. Between 20 and 30 km, occasionally large monthly 
variations in location and thickness of these layers are observed. The 
data suggest the presence of a quasi -permanent aerosol scattering layer 
of 5 km thickness centered at 20 km. 
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ALTITUDE 



AEROSOL BACKSCATTERING COEFFICIENT per meter per sterodlan 


FIGURE I TIME SERIES OF AEROSOL BACKSCATTERING COEFFICIENT f A , VERSUS HEIGHT COMPUTED FROM 
LIDAR BACKSCATTER OBSERVATIONS MADE AT MENLO PARK, CALIFORNIA 

Error bars represent ± one standard deviation of the data sample. Dashed line represents constant mixing ratio. 

(Drawn at 1 percent of the molecular; backscattering coefficient) 



COMPARISON OF LIDAR AEROSOL MEASUREMENTS 


AND DIRECT PARTICULATE SAMPLING 

Frederick G. Fernald, Burton G. Schuster, 
and Charles L. Frush • 

National Center for Atmospheric Research* 

Boulder, Colorado 

ABSTRACT 

Lidar data were collected at Boulder , Colorado , on two evenings ; 308 
shots between 21x3k and 00:33 (MST) on the night of September 12-13, and 152 
shots between 21:U3 and 22:51 (MST) on September 15. The averaged lidar re- 
turns were normalized with respect to the real molecular atmosphere as deter- 
mined by Denver, Colorado, radiosondes. Error limits of ± 0.67c become dis- 
tinguishable from the normalized returns only above 25 km. 

The data for September 12 were normalized to the real atmosphere by de- 
termining the level at which the lidar return (average of all 308 shots) was 
a minimum compared to the molecular atmosphere and setting the ratio to 1.0, 
that is by assuming that there were no aerosols present in the atmosphere at 
that level. For the 12th of September this level was 9*5 km. The same pro- 
cedures, though, could not be successfully applied to the September 15 data 
for it appears that significant amounts of aerosols were present at all 
heights (between 7*5 and 22 km) to which this normalization method was 
applied. The data for the 15th were , therefore , adjusted so that they agreed 
with the results of the 12th of September in the l6 to 2k km height interval. 

The corresponding aerosol backscattering cross sections have been plotted 
superimposed upon the University of Wyoming balloonsonde aerosol counter data. 
The agreement among these three sets of data in the stratosphere (above lL km) 

* The National Center for Atmospheric Research is sponsored by the National 
Science Foundation. 
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is very good considering their spatial and temporal differences. A cirrus 
cloud layer appearing between 10 and 12 km on the 15th of September has not 
been plotted 5 since the relatively large backscattering cross sections of 
these sub-visible cirrus particles greatly exceeded the scale of the graph. 

Perhaps most important of all, it appears that an integration from the 
tropopause to 2k km would produce results within a few percent for the three 
sets of data. 
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GLOBAL RECONNAISSANCE OF STRATOSPHERIC AEROSOL 
BY AIRBORNE LIDAR 

Burton G. Schuster , Frederick G, Fernald, 
and Charles L. Frush 

National Center for Atmospheric Research* 

Boulder , Colorado 

ABSTRACT 

As part of the Climatic Impact Assessment Program (CIAP), a series of 
airborne lidar measurements have been performed to eventually determine , on a 
global scale, the extent and degree of particulates in the stratosphere. Thus 
far, two missions have been flown. The first, starting 25 February 1973, went 
from Albuquerque, New Mexico - Goose Bay, Labrador - Fairbanks, Alaska - 85°N 
latitude - Hawaii - equator. The second mission, starting 1 April 1973, went 
from Albuquerque, New Mexico - Goose Bay, Labrador - Iceland - Madrid, Spain - 
Albuquerque. Several local Iceland missions were flown to determine if the 
Helgafell volcanic activity had penetrated into the stratosphere. 

The items to be discussed will include a general discussion of the overall 
global mission, the experimental interface of the lidar and the Air Force 
NC-135 aircraft, and a discussion of the preliminary results to date. These 
latter include nacreous clouds and associated low temperatures in mountainous 
northern latitudes; and layers up to 30 km in the equatorial regions. 


*The National Center for Atmospheric Research is sponsored by the National 
Science Foundation. 



COMPARISON OF LIDAR AND IN-SITU 


MEASUREMENTS OF STRATOSPHERIC AEROSOLS 

by 

S. H. Melfi, G. B. Northam, M. P. McCormick (NASA Langley Research Center) 

and 

J. M. Rosen, T. J. Pepin-, D. H, Hofmann (University of Wyoming) 

ABSTRACT 

This paper will present the results of a comparative study conducted in 
Laramie, Wyoming, during the summer and fall of 1972, as part of the Depart- 
ment of Transportation's Climatic Impact Assessment Program (CIAP). The 
study included independent, and nearly simultaneous, measurements of strato- 
spheric aerosols using a LIDAR system and a balloon-borne in-situ particle 
counter. 

The LIDAR provides a remote measurement of volume backscatter (aerosols 
and molecules) in a narrow wavelength region centered at the ruby wavelength 
( 69 ^- 3 ^); whereas the balloon-borne in-situ counter measures aerosol concentra- 
tion by counting aerosols greater than - 0.30 ]im in a diameter as they are 
pumped through a chamber and scatter white light forward into photo-detec- 
tors . 

The comparison of measurements that will be discussed using the two 
techniques involves formulating the LIDAR data so that it is compatible with 
the counter data. The formulation includes separation of the scattering due 
to aerosols from the total and displaying this in terms of aerosol scattering 
function. Aerosol scattering function is proportional to aerosol concentra- 
tion if the aerosol parameters, such as size distribution and composition, are 
constant with altitude. 

In separating the aerosol scattering from the total, the need for real 
atmospheric number density over the Standard Atmosphere is also discussed. 
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STRATOSPHERIC AEROSOL LIDAR MEASUREMENTS 


OVER JERUSALEM 

A. Cohen and M. Grafter 
Department of Atmospheric Sciences 
The Hebrew University of Jerusalem, Israel 

ABSTRACT 

The laser radar system , consisting of a ruby laser as the emitter, has 
been used to map the aerosol layer centered at 20 km. Scattering measure- 
ments were made in two polarizations and changes as a function of time in the 
number density as well as in the size distribution were detected. 

The size distribution variations were observed by the analysis of the de- 
polarization of the backscattered light. 

The results were compared to those obtained in different sites over the 
whole globe and the height of the maximum concentrations vs. latitude and 
season were plotted. 

Regular measurements of the stratospheric aerosol profile are being car- 
ried out in order to detect future concentration variations due to volcanic 
erruptions or other possible dust sources. The profile describes the altitude 
dependence of the aerosol concentrations between 10 - 30 km, by use of a 
multichannel analyzer. 
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LASER RADAR OBSERVATIONS ON THE MESOSPHERE AND LOWER THERMOSPHERE 

OVER KINGSTON 5 JAMAICA 


Ly 

G, S. Kent and W. Keenliside 
Department of Physics, University of the West Indies 

ABSTRACT 

The Mark II laser radar system has been in use on a routine basis for 
nearly three years. Typically three continuous operational periods of 8-10 
hours are made each month., and the variations of the signal scattered at 
heights from 60-100 kilometres studies. The normal integration period used is 
approximately 30 minutes and the relative accuracy of measurement of the re- 
ceived signal;, under good atmospheric conditions, varies from about 1 % at 70 
km to 15 % at 100 km. The following properties of the atmosphere in this 
height range have been established. 

1. (a) The solar diurnal tidal oscillation has a dominant wavelength of 15 

km, which is within the range of wavelengths predicted but which 
cannot be assigned to any specific mode. The wave is strongly at- 
tenuated with increasing altitude. Seasonal variations are observed 
in wavelength and phase; spring and autumn showing the most regular 
behaviour . 

(b) The solar semi-diurnal tide is found to agree well with theoretical 
predictions. It does not appear to be greatly attenuated with alti- 
tude and is stronger than the diurnal tide above about 90 km. Below 
this height the diurnal tide is dominant. Seasonal changes in phase 
are observed. 

2. A seasonal variation in atmospheric density is observed. Maximum 
density is observed in winter and minimum in summer; the annual 
range is approximately 10 % at a height of 90 km. 
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There is no regular mesospheric dust or aerosol layer present 
over Jamaica. 

Considerable day-to-day variability in behaviour is observed, 
both in the background density and the tidal oscillations. 

Good agreement has been obtained between the characteristics of 
the atmospheric tides as observed by the laser radar system and by 
meteor radar and short-wave radio sounding at the same location. 



MEASUREMENT OF RESONANCE SCATTERING CROSS SECTION OF SODIUM D LINES 
AND LASER RADAR DETECTION OF SODIUM LAYER IN THE UPPER ATMOSPHERE 

BY A TUNABLE DYE LASER 

M. Jyumonji, T. Kobayasi, H. Inaba, T. Aruga* and H. Kamiyama# 

Research Institute of Electrical Communication, Tohoku University 

Sendai, Japan 

^Institute of Geophysics, Faculty of Science, Tohoku University 

Sendai, Japan 
ABSTRACT 

The recent development of tunable dye laser technology is promoting a new 
approaches to high resolution spectroscopy, atmospheric research and air pollu- 
tion analysis. This paper describes some output characteristics of a tunable 
dye laser and its application to the measurement of resonance scattering cross 
section. Furthermore, the measurement of the sodium layer in the upper 
atmosphere with a resonance scattering laser-radar system performed by our 
laboratory is reported. 

A flashlamp-pumped dye laser was designed and constructed for the resonance 
scattering study. A combination of a linear flashlamp and an elliptical cylin- 
drical reflector was employed to pump uniformly the Rhodamine 6G solution for 
obtaining high frequency stability. The flashlamp is driven by a 0.5 yF low 
inductance capacitor, which produced a pumping light with a fast rise time 
shorter than 0.3 yS. Spectral narrowing of the oscillation width and tuning 
of the wavelength are accomplished with a tilted Fabry-Perot filter mounted 

between a diffraction grating and a mirror. The laser output is typically about 

° . o 

5 mJ and the spectral width varies from 0.01 A to 0.04 A as the flashlamp input 
increases. 

The resonance scattering cross section of sodium D lines were measured 
carefully by irradiating the dye laser beam in a 3-cm diam spherical sodium 
vapor cell. Table 1 shows the theoretical and experimental values of the 
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Table 1. 


Na I line 

D 1 

D 2 

Measured direction 

0,7T 

tt/2 

0,TT 

tt/2 

(da/dh) 

do -13 

cm^/sr ) 

Experimental 

- 

3-3 

- 

6 .H 

Theoretical 

3*1 

3.1 

6.5* 

7 . 7 ** 

6.9*** 

5 . 9 * 

5 .H** 

5 . 7 *** 


Calculated on the basis of the two-photon 
process including the hyper-fine structure. 
**Calculated on the basis of the two-photon 
process including the fine structure. 
***Calculated from the absorption crosssection. 


differentical cross sec- 
tion for the resonance 
scattering of the sodium 
and Dg lines. The 
theoretical values were 
calculated on the basis of 
two different interaction 
schemes; l) absorption, 
and 2) two-photon process. 

The cross sections measured 
at a right angle to the laser 
_13 and 


beam were 3*3X10 
6.Hxl0 -13 cm^/sr for 
and Dg lines, respectively. These values indicate a good agreement with 
theoretical values. 

The laser radar system using this dye laser as a transmitter was 
installed at the Zao Upper Atmospheric Observatory, Tohoku University 
(lU0.5 O E, 38.3°N) in order to detect the upper atmospheric sodium layer. A 
50-cm diam Cassegrain telescope and a 10-channel photon counter were used for 
receiving the resonance backscatter. Seasonal variation of the structure of 
the sodium layer has been observed during 1972. Figure 1 shows one of the 
typical observed results performed on June lH-15, 1972, indicating the loca- 
tion of the sodium layer. As an interesting result of these measurements. 


Table 2. 


Date 

Time 

Number density (l 03 cm - 3 ) 

— 

Column 

Altitude (km) 


density 



Peak 1 

Peak 2 

Peak 3 

(l 09 /cm 2 ) 

Jun. lH 

2220-0010 

ITU 

0.2 



(90-95) 

(120-125) 


X . X 

Oct. 7 

0030 - 0 H 35 

1.0 

1.2 

0.9 

3.1 

(90-95) 

(100-105) 

(120-125) 

Oct. 8 

2050-2200 

1.8 

2.3 

1.2 


( 75 - 80 ) 

(90-95) 

(120-125) 

3 * 1 

Oct. 10 

O 

O 

1 

O 

-3" 

CM 

O 

2.9 

1.6 

1.7 

7-2 



( 90 - 95 ) 

(uo-115) 

(120-125) 

Oct. 13 

03H0-0H30 

0 .H 

0.5 


l.H 

( 90 - 95 ) 

(115-120) 



Ho 





for instance, an increase in the N abundance was evidently observed just 

SL 

after the Jacobini meteoric shower on Oct. 8, 1972. 

Table 2 illustrates the change of sodium density distribution before and 
after this meteoric shower. This result may give one of possible explanations 
concerning the origin of the atmospheric sodium. 

Other experimental results will be reported in detail. 

1) T. ¥. Hansh et al, Phy. Rev. Letters, 27 , 707 (1971). 

2) M. R. Bowman et al, Proc. Radio and Electronic Engineer, 39 , 29 (1970). 

3) M. Jyumonji, T. Kobaysi, H. Inaba, T. Aruga, and H. Kamiyama, Proc. 9th. 
Symp. on Space Tech, and Sci. Tokyo, 123 (1970)* 
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FLUORESCENCE OF AMBIENT ATMOSPHERIC 
AEROSOLS AND ITS IMPLICATIONS FOR REMOTE 
DETECTION OF GASEOUS ATMOSPHERIC POLLUTANTS 

Jerry Gelbwachs and Milton Birnbaum 
The Aerospace Corporation 
Electronics Research Laboratory 
Los Angeles , California 90009 

ABSTRACT 

The visible fluorescence of aerosols in the ambient atmosphere has been 

1 

monitored m situ using ew argon ion laser excitation. Fluorescence was ob- 
served at 90° in bands of 500 and 1000 X over the spectral range of O .56 to 
0.8l ym. Aerosol and nitrogen dioxide (NO ) fluorescence and Mie scattering 
were monitored in the ambient air in practically real time. The buildup and 
decline of the aerosol levels and N0^ concentrations were similar (Fig. l). 

The observed aerosol fluorescence can constitute a limitation upon the 

ability of a LIDAR to measure pollutant concentrations. The aerosol data in 

conjunction with measured Raman scattering from nitrogen permits estimates as 

2 

to the extent of aerosol interference to remote Raman measurements. For ex- 
ample, under severe aerosol loading in the Los Angeles basin, the fluorescence 
of aerosols generates a signal equivalent to 600 ppm of SO^ when detecting the 
Stokes wavelength (bandwidth of 35 X) excited at L880 

This limitation can be largely overcome by utilizing a two wavelength 
excitation method where the aerosol fluorescence is approximately constant 
while the molecular signals differ. For many LIDAR applications the back- 
ground radiation exceed the molecular signals. Under this condition it is 
shown that the degradation in signal to noise ratio with aerosol fluorescence 
compared to the identical situation without aerosol fluorescence is directly 
proportional to (l-8); where 8 is the ratio of the molecular signals at the 
two wavelengths, 8 less than unity. 

We have demonstrated the effectiveness of the two wavelength method by 
in situ measurements of atmospheric NO^ levels in the presence of ambient 
aerosol fluorescence. Excitation was provided at 4965 2. and 5017 2. The 
fluorescence signals from NO^ monitored at O.65 to 0.75 ym were twice as 

42 



large when excited H 965 £ as compared to 5017 £ excitation while the signals 
from the aerosol remained constant. The two wavelength scheme accurately mea^ 
sured the ambient levels of several parts per hundred million (pphm). 
Therefore, by application of appropriate techniques, low pollutant concentra- 
tions may still be determined by LIDAR in the presence of aerosol fluores- 
cence . 


REFERENCES 

1. Apparatus similar to that of M. Birnbaum, J. A. Gelbwachs , A. W. Tucker, 
and C. L. Fincher, "Fluorescence Determination of Atmospheric NOg," 
scheduled for publication in OPTO-ELECTRONICS , May, 1972. 

2. D. G. Fouche and R. K. Chang, "Relative Raman Cross Section for Ng, 0' 2 , 

’ CO, C0 2 , S0 2 and H^," APPL. PHYS. LETTERS, Vol. 18 , pp 579-580, 1971. 
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S0 2 SPECTROSCOPY WITH A TUNABLE UV LASER 

W. W. Morey, C-. M. Penney and M. Lapp 
Optical Physics Branch 

General Electric Co. , Corporate Research and Development 

Schenectady, New York 

ABSTRACT 

A portion of the fluorescence spectrum of SO^ has heen studied using a 
narrow wavelength doubled dye laser as the exciting source. One purpose of 
this study is to evaluate the use of SO^ resonance re-emission as a probe of 
SO^ in the atmosphere. 

When the SO^ is excited by light at 300.2 nm, for example, a strong re- 
emission peak is observed which is Stokes-shifted from the incident light 
wavelength by the usual Raman shift (the symmetric vibration frequency 
1150.5 cm -1 ). 

The intensity of this peak is sensitive to small changes (.01 nm) in the 
incident wavelength. Measurements of the N^ quenching and self quenching of 
this re-emission have been obtained. Preliminary analysis of this data indi- 
cates that the quenching is weak but not negligible. 

The dye laser in our system is pumped by a pulsed N^ laser. Tuning and 
spectral narrowing are accomplished using a telescope-echelle grating combina- 
tion. In a high power configuration the resulting pulses have a spectral 
width of about 5 x 10 nm and a time duration of about 6 nsec. The echelle 
grating is rotated by a digital stepping motor, such that each step shifts the 
wavelength by 6 x 10 ^ nm. 

In addition to the tunable, narrow wavelength uv source and spectral 
analysis of the consequent re-emission, the system also provides time resolu- 
tion of the re-emitted light to 6 nsec resolution. This capability is being 
used to study the lifetime of low pressure SO^ fluorescence at different wave- 
lengths and pressures. 

* Supported by NASA contract NAS1-11624 
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INELASTIC LIGHT SCATTERING PROCESSES 



Daniel G. Fouche* and Richard K. Chang 

/ 

Department of Engineering and Applied Science 
Yale University 

New Haven, Connecticut 06520 
ABSTRACT 

Five different inelastic light scattering processes will he denoted by, 
ordinary Raman scattering (ORS), resonance Raman scattering (RRS), off-reso- 
nance fluorescence (ORF), resonance fluorescence (RF), and hroad fluorescence 
(BF), A distinction between fluorescence (including ORF and RF) and Raman 
scattering (including ORS and RRS) will he made in terms of the number of in- 
termediate molecular states which contribute significantly to the scattered 
amplitude , and not in terms of excited state lifetimes or virtual versus real 
processes . 

The theory of these processes will be reviewed, including the effects of 
pressure, laser wavelength, and laser spectral distribution on the scattered 
intensity. 

The application of these processes to the remote sensing of atmospheric 
pollutants will be discussed briefly. It will be pointed out that the poor 
sensitivity of the ORS technique cannot be increased by going toward resonance 
without also compromising the advantages it has over the RF technique. 

Experimental results on inelastic light scattering from I vapor will be 
presented. As a single longitudinal mode 51^5 A argon-ion laser line was 
tuned away from an I absorption line, the scattering was observed to change 
from RF to ORF. The basis; of the distinction is the different pressure de- 
pendence of the scattered intensity. Nearly three orders of magnitude en- 
hancement of the scattered intensity was measured in going from ORF to RF. 
Forty-seven overtones were observed and their relative intensities measured. 

The ORF cross section of I compared to the ORS cross section of N. was found 
6 2 2 
to be 3 x 10 , with I^ at its room temperature vapor pressure. 

* Current address: M.I.T. Lincoln Laboratory, Lexington, MA 02173. 


46 



ATMOSPHERIC TEMPERATURE FROM RAMAN SCATTERING 


M. Lapp and C. M. Penney 
Optical Physics Branch 

General Electric Co., Corporate Research and Development 
Schenectady, New York 

ABSTRACT 

Raman scattering signatures are functions of the scattering-gas tempera- 
ture, and are therefore of potential use for practical atmospheric temperature 
probes. The method described uses either ratios of pure rotational Raman 
scattering intensities for air utilizing various spectral bandpasses, or, al- 
ternatively, ratios of air-rotational to nitrogen-vibrational scattering in- 
tensities. Three aspects of work relating to the development of such probes 

1 

are discussed in this presentation: (l) fundamental absolute Raman data, 

(2) construction of air spectra from these data, and (.3) temperature-sensi- 

2 

tivity of the signature. 

The fundamental data described are the absolute rotational and vibration- 
al scattering cross sections. Recent measurements in this laboratory of rota- 
tional cross sections for , 0^ , and C0^ are emphasized, as are their use in 
predicting absolute magnitudes of Raman scattering signals. Next is described 
the computation of air rotational Raman spectra as a function of temperature 
calculated through use of the experimentally-measured cross sections. The 
spectra are based additively upon nitrogen and oxygen contributions, since 
pure rotational Raman scattering from water vapor is very weak. 

Finally, the sensitivity of the scattering intensity ratios to tempera- 
ture is explored as a function of choice of spectral bandpass for the monitored 
rotational Raman scattering. Various compromises will be discussed which must 
be made in choosing bandpasses appropriate for specific purposes and experi- 
mental conditions. 
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This work was supported in part by NASA Lewis Research Center on Contract 
NAS3-15825 . A more complete description of this material can he found in 
"Absolute Intensity and Polarization of Rotational Raman Scattering from 
Ng, 0^ and COg" by C. M. Penney, R. L. St. Peters, and M. Lapp, MSA Re- 
port No. 0R-121091, 2 January 1973. 

2. This work was supported in part by Aerospace Research Laboratories, 
Wright-Patterson Air Force Base on Contract F33615-71-C-1867 . A more 
complete description of this material can be found in "Application of 
Light-Scattering Techniques for Measurements of Density, Temperature, and 
Velocity in Gas-dynamics," by M. Lapp, C. M. Penney, and J. A. Asher, 
Aerospace Research Laboratories, WPAFB, Report No. ARL 73-004-5 (in press). 
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REMOTE MEASUREMENT OF ATMOSPHERIC 
TEMPERATURES BY RAMAN LIDAR 


by 

Jack A. Salzman and Thom A. Coney 
National Aeronautics and Space Administration 
Lewis Research Center, Cleveland, Ohio 

ABSTRACT 

The Raman shifted return of a lidar, or optical radar, system has been 
utilized to make atmospheric temperature measurements. These measurements 
were made along a horizontal path at temperatures between -20°C and +30°C 
and at ranges of about 100 meters. The temperature data were acquired by re- 
cording the intensity ratio of two portions of the Raman spectrum which were 
simultaneously sampled from a preset range. 

The lidar unit employed in this testing consisted of a It joule-lOppm laser 
operating at 694.3 nm, a 10-inch Schmidt-Cassegrain telescope, and a system 
of time-gated detection and signal processing electronics. The detection 
system processed three return signal wavelength intervals - two intervals along 
the rotational Raman scattered spectrum and one interval centered at the 
Rayleigh-Mie scattered wavelength. The wavelength intervals were resolved by 
using a pellicle beam splitter and three optical interference filters. Raman 
return samples were taken from one discrete range segment during each test shot 
and the signal intensities were displayed in digital format. The Rayleigh-Mie 
return was monitored continuously through standard oscilloscope display tech- 
niques . 

The test site utilized to evaluate this measurement technique encompassed 
a total path length of 200 meters. Major components of the test site included 
a trailer-van housing the lidar unit, a controlled environment test zone, and a 
beam terminator. The control zone which was located about 100 meters from the 
trailer was 12 meters in length, 2.4 meters in diameter, and was equipped with 
hinged doors at each end. The temperature of the air inside the zone could be 
either raised or lowered with respect to ambient air through the use of 



infrared heaters or a liquid-nitrogen cooling system. Conditions inside the 
zone were continuously monitored with a thermocouple rake assembly. The test 
path length was terminated by a 1.2 meter square array of energy absorbing 
cones and a flat black screen. 

Tests were initially conducted at strictly ambient conditions utilizing 
the normal outside air temperatures as a test parameter. These tests provided 
a calibration of the Raman intensity ratio as a function of temperature for 
the particular optical-filter arrangement used in this system while also pro- 
viding a test of the theoretical prediction formulated in the design of the 
system. Later tests utilized zone temperatures above and below ambient to 
provide temperature gradient data. These tests indicate that ten shots, or one 
minute of data acquisition, from a 100 meter range can provide absolute temper- 
ature measurements with an accuracy of + 3°C and a range resolution of about 
5 meters. Because this measurement accuracy compares well with that predicted 
for this particular unit, it is suggested that a field-application system could 
be built with significant improvements in both absolute accuracy and range. 
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THE USE OE TUNABLE LASERS IN INFRARED 
ATMOSPHERIC SPECTROSCOPY* 


P. L. Kelley 

Lincoln Laboratory, Massachusetts Institute of Technology 
Lexington, Massachusetts '02173 

ABSTRACT 

A review will be given of tunable sources of coherent infrared radiation - 
with particular emphasis on lead salt diode lasers. Characteristics such as 
linewidth, power output and mode tuning will be discussed. 

High resolution tunable laser studies have been carried out for molecules of 
interest in atmospheric studies, such as H^O , SOg , NO, CO and C 2 H^. Measure- 
ments have been made on widths, strengths, shapes and relative positions of 
absorption lines. Dicke narrowing, hyperfine structure, A-doubling, and Stark _ 
and Zeeman tuning have been studied. 

Air pollution measurements can be made with tunable infrared lasers ; these 
include point sampling, long path atmospheric transmission, and heterodyne 
spectroscopy of remote emission sources. Online automobile exhaust measure- 
ments of ethylene, in situ smokestack measurements of SOg , and long path 
atmospheric measurements of ethylene will be described. 


*This work was sponsored by the Department of the Air Force. 



COMPARISON OF LIDAR METHODS FOR REMOTE MEASUREMENT 
OF AIR POLLUTANTS 

4 * 

M. L. Wright, E. K. Proctor and E. M. Liston” 

ABSTRACT 

This paper presents quantitative comparisons of several single-ended lidar 
techniq^les for the remote measurement of gaseous pollutants. These techniques 
are divided into two groups. The first group is Based on the measurement of 
energy scattered directly by the gas of interest. The gaseous scattering 
processess considered are ordinary fluorescence, resonance fluorescence (also 
called resonance scattering), Raman scattering, and resonant (or nearly 
resonant) Raman scattering. The second group is based on the measurement of a 
characteristic differential absorption produced by the gas of interest at two 
discrete wavelengths, using energy scattered back toward the receiver by a 
remote reflector other than the gas of interest. The remote reflector may be 
intermixed with the gas of interest, as is the case with aerosols and atmospher- 
ic gases (principally nitrogen) , or they may be fixed reflectors such as ter- 
restrial objects or retroflectors . 

The detectability of a given material will depend on the magnitude and 
characteristics of the optical interaction with that material. The main 
characteristics of interest are the cross section, the response time, and the 
spectral response of the material relative to both the transmit and receive 
functions of the lidar. These characteristics and their implications for re- 
mote sensing will be reviewed for the four direct scatter processes and for 
the differential absorption technique. 

The characteristic behavior of the direct backscatter technique is dif- 
ferent from the differential absorption technique with respect to sensitivity, 
concentration of material, and the effect of range. For these reasons, the 
direct backscatter processes cannot be compared directly to the differential 
absorption technique. The two techniques can be compared for specific material 
and system configurations, however. This paper describes specific lidar system 
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configurations and gives the calculated performance level for these systems in 
both the direct backscatter and differential absorption modes for a wide variety 
of pollutant monitoring situations. 

The results of this comparison of techniques indicate that the differen- 
tial-absorption lidar technique can provide adequate range and sensitivity for 
a wide variety of pollution monitoring applications involving a number of in- 
teresting pollutant materials. No other single technique appears to provide 
these capabilities for such a wide range of materials. 


^Systems Techniques Laboratory 
tRadio Physics Laboratory 
^Atmospheric Sciences Laboratory 
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AN ERROR ANALYSIS OF THE DIFFERENTIAL ABSORPTION TECHNIQUE 


Richard M. Schotland Met 0 19 
New York University* 

ABSTRACT 


The measurement of the spatial distribution of atmospheric gases by means 
of a lidar operating in the absorption mode is presently being investigated by 
a number of groups. In this technique a pulsed laser is tuned to the center of 
an isolated absorption line of the gas under consideration and a measurement is 
made of the backscattered radiation as a function of distance. A similar 
measurement is then taken with the laser tuned to an adjacent window region so 
that the absorption and scattering processes of the atmosphere can be separated 
from those of the gas. The spatial distribution of the gas is obtained by 
differencing the logarithm of the ratio of the two measurements of back scat- 
tered power in such a manner as to obtain the differential absorption of the 
laser radiation by the gas with respect to the lidar range. This technique 
has been referred to by the mnemonic DASE (Differential Absorption of Scat- 
tered Energy). The prupose of this paper is to present an error analysis of 
the DASE technique with particular reference to the measurement of atmospheric 
water vapour . 

The power captured by a lidar receiver is given by the expression: 


' rV 


cJ v A 

2R 2 


B,,exp - 2 /f (Y v + 2p K v ) dr 


V 


V' 


(1) 


An expression for p averaged over a distance AR can be obtained if equation 
(l) is evaluated at ranges R and R + AR for back scatter measurements made at 
frequencies corresponding to line center (v=L) and to the window (v=w). 


P (R) 


l 

2KAR 
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P rL (R) P 'rW (R) 

1 , TLi - I rW + rp# + -d* 

n p (R+AR) n p T tt (R+AR) 1 B 
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where 


T* = -2 



. 6 l (R+AR) B l (R) 

B * * 111 b’ w (h+ab) - ln sys) 


The prime notation indicates that the window measurement is made at a time dif- 
ferent than that of the line measurement . 

p (R) can he determined from (2) if T* and B* are known. It can he shown that 

y T ~ y TT and $ T ~ 3 TT for frequency shifts of L - W. Further, it can he demon- 
L W L W 2 

strated that if the measurement pairs are made within 10 seconds then 
Y x ~ y'tt and £ T ~ 3*. Consequently B* and T* can he taken to he zero for a 
suitably designed lidar system. 

The variance in the deduced quantity p due to the nature of the measure- 
ment process can he shown to he approximately equal to the expression: 



p 2 k 2 ar 2 

V 



(3) 


The first "bracketed term represents that portion of the fractional variance of 
p which is accounted for by uncertainties in the knowledge of K^. The un- 
certainty in arise from unknown variations of temperature and pressure along 
the path of the lidar path as well as in the knowledge of V. 

_ 2 




AM 2 

(2 nta ) 2 + AM 2 

Jj 


( 4 ) 


Here a and b are of the order of unity. AM gives the uncertainty in 

the axial mode number for a laser whose cavity length is nt , tuned to a Lorentz 

line whose half width is a T . 

Jb 

Substitution of typical parameters into equation (k) leads to an estimate 
of the fractional standard deviation of to be approximately 
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The second "bracketed term of equation (3) involves which is the ’noise 

rr 

to signal ratio of the measured "back scattered power. The information content 
of the measurement is equal to the absorption of the laser radiation over the 
path AR. Therefore the uncertainty in the deduced value of p will depend upon 
the noise to signal value of the measurement divided by the power loss in AR. 

AR is normally taken so that the product p K AR is approximately 0.1. There- 
fore the fractional standard deviation of p is of the order of 10 times the 
noise to signal ratio of the power measurement. 

An application of the analysis presented will be made to the measurement 
of tropospheric water vapor by the DASE technique. 


#0n Sabbatical leave at Meteorological Office 4 Bracknell 4 U. K. 



LASER RADAR STUDY USING RESONANCE ABSORPTION 
FOR REMOTE DETECTION OF AIR POLLUTANTS 


Takashi Igarashi 

Radio Research Laboratories , Ministry of Posts 
and Telecommunications , Koganei , Tokyo , Japan 

ABSTRACT 

A laser radar using' resonance absorption has an advantage of increased 
detection range and sensitivity compared with that achieved by Raman or 
resonance back scattering . 

In this paper 9 new laser radar system using resonance absorption is 

proposed and results obtained from this laser radar system are discussed, 

. o o 

N0 2 , SO^ gas has an absorption spectrum at 4500 A and 3000 A 
respectively as shown in Fig, 1, A laser light including at least a set 
of an absorption peak A^ and a valley A^ is emitted into a pollutant 
atmosphere. The light reflected with a topographical reflector or an 
atmospheric Mie scattering as distributed reflectors is received and 
divided into two wavelength components A , A^. From the ratio of these two 
components Y (x) , the distribution of the pollutant gas n(x) can be calculated 

■by 


n (x) - i SJLfei 

2 (c^ - a 2 ). d x 

5 where 0^5 correspond to the absorption cross sections of the pollutant 
gas at A^ and A 2 respectively. 

The laser radar system used in the investigation is shown in Fig,. 2 

and consists of a dye laser transmitter, an optical receiver with a special 

monochrometer and a digital processer. Table 1 shows the molecular constants 

of NO , and SO and the dye laser used in this experiment. 

2 d 
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o 

no 

S° 2 

A 1 

(A) 

1 + 1+82 

3001 

*2 

1+1+69 

3011+ 

2 

a i “ °2 (°m ) 

2.0 x 10" 19 

2.5 x 10~ l8 

\ - a 2 (A) 

13 

13 

Dye Material 

calcein blue 

1 

SHGed 

Ehodamine B 

0 

Spectrum width (A) 

30 

30 

Pumping Source 

Sum of SHGed 
and 

fundamental 
Nd: YAG 

SHGed 
Nd: YAG 


Table 1 Molecular constants of NO^ and SO^ and dye laser used 

In this system, the absolute concentration of the pollutant gas can 
be measured in comparison with a standard gas cell. The concentration of 
UO^, SO^ as low as 0.1 ppm have been measured at 100 m depth resolution. 

For a 1 mJ laser output , the observable range of this system achieved up 
to 300 m using the distributed Mie reflector. 

The capability and technical limitation of the system will be discussed 
in detail. 
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Fig.l ‘Absorption spectrum of. NC^ and SO2 



Fig. 2 Schematic diagram of laser radar system 
using resonance absorption for detection of air 
pollutants. 
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DETERMINATION OF SPATIAL DISTRIBUTION 
OF AIR POLLUTION BY DYE LASER MEASUREMENT OF 
DIFFERENTIAL ABSORPTION OF ELASTIC BACKSCATTER 


S.A. Ahmed and J.S. Gergely 
Electrical Engineering Department 
The City College of New York 
New York , New York 10031 
ABSTRACT 

This paper presents the results of an analytical study of a lidar system 
which uses tunable organic dye lasers to accurately determine spatial distri- 
bution of molecular air pollutants. Also described will be experimental 
work to date on simultaneous multiwavelength output dye laser sources for 
this system. 

Basically the scheme determines the concentration of air pollutants 
by measuring the differential absorption of an (at least) two wavelength 
lidar signal elastically backscattered by the atmosphere. Only relative 
measurements of the backscattered intensity at each of the two wavelengths, 
one on and one off the resonance absorption of the pollutant in question, 
are required. 

The various parameters of the scheme are examined and the component 
elements required for a system of this type discussed, with emphasis on the 
dye laser source. Potential advantages of simultaneous multiwavelength 
outputs are described. The use of correlation spectroscopy in this context 
is examined. Comparisons are also made for the use of infrared probing 
wavelengths and sources instead of dye lasers . 

Estimates of the sensitivity and accuracy of a practical dye laser 
system of this type, made for specific pollutants, show it to have inherent 
advantages over other schemes for determining pollutant spatial distribution. 
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REMOTE SENSING OF METHANE USING AN ERBIUM/YAG LASER: 

A FEASIBILITY STUDY 

by 

Kenneth 0* White 
Gerald T. Wade 

Atmospheric Sciences Laboratory 
US Army Electronics Command 
White Sands Missile Range, New Mexico 

Stuart A,. Schleusener 
Electrical Engineering Department 
New Mexico State University 
Las Cruces, New Mexico 

ABSTRACT 

A program is underway to determine the feasibility of a remote sensing 
scheme for the detection of methane based on a coincidence "between the emission 
of the Erhium/YAG laser and a methane absorption line. This recently discovered 
coincidence occurs with the R(6) line in the 2v^ overtone band of methane. The 
methane line is centered at 16^5. lnm and is about O.lnm wide; the Erbium/YAG 
laser emission is centered at l61|l+.9nm and is about 0.6nm wide. Theory indi- 
cates that the absorption line is composed of four or more components, and 
initial measurements indicated that there was structure in the line. 

Two remote sensing techniques are being considered: the differential 

absorption transmission method and the differential absorption backscatter 
technique. In order to obtain results, the line parameters for the methane 
absorption line are needed. However, very little information is available on 
this absorption line. The research program consists of an experimental part, 
to determine the spectral line parameters, and a theoretical part, to deter- 
mine performance* The experimental program consists of two parts: a measure- 

ment of the absorption coefficient for the line, and measurements to determine 
the structure of the line and the positions , halfwidths , and strengths of the 
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components. The theoretical program was also divided into two sections. The 
first was a study using the lest available data. The second part uses more 
accurate experimental data as it becomes available. The results of the study 
using the best available data indicate for the differential absorption trans- 
mission method that if the receiver system can detect a 1 % change in trans- 
mission, then a 5% change in the methane concentration can be detected with a 
range of 100m to 8km, depending on the value of the absorption coefficient. 
Results of this study have been plotted as change in transmission versus range 
for various changes in the methane concentration for several values of the 
absorption coefficient. The bounds on the parameters were: minimum detectable 

change in transmission, 0.01 to 100%; range, 1m to 10km; change in methane 
concentration , 1 to 5000 %, and absorption coefficient, 0.025 to I.65 ppm \:m \ 

For the differential absorption backscatter technique, preliminary work 
consisted of comparing the available methane line parameters with those of the 
0^ and H^0 lines in the region of the ruby laser emission. A detection scheme 
based on the ruby laser is feasible for the detection of 0^ and H^O. The data 
indicates that the absorption coefficient for the absorption lines of methane 
are up to an order of magnitude larger. However, the wavelength dependence 
of the scattering and the detector parameters available at 1.6 micrometers 
combine to decrease the available signal by 2 to 3 orders of magnitude. 

The above results indicate that the measurement of methane using an 
Erbium/YAG laser and the differential absorption transmission technique is 
feasible. However, the use of the differential absorption backscatter method 
appears marginal unless sophisticated detection methods, such as photon count- 
ing, are used. Results of the experimental program and the second part of the 
theoretical program will also be presented. 
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ON ATMOSPHERIC TEMPERATURE MEASUREMENT 
BY LIDAR DIFFERENTIAL ABSORPTION 

by 

James B. Mason 

Atmospheric Sciences Laboratore 
US Army Electronics Command 
White Sands Missile Range, New Mexico 

ABSTRACT 

The differential absorption of laser energy in the atmosphere by differ- 
ent lines within a molecular rotational band can be measured to obtain the 
relative strength of the lines. Since line strengths are directly related to 
local temperature via the partition function, a method for measuring tempera- 
ture is suggested. 

By considering the propagation of lidar energy at both absorbed and non- 
absorbed wavelengths , a simple relationship has been derived via the lidar 
equation which expresses the temperature of the absorbing gas in terms of the 
backscattered signals only. It is found that the temperature as related to 
the ratios of the scattered powers and the spatial derivatives of these ratios. 

A computer program incorporating resonance absorption has been used to 
calculate the lidar returns from model atmospheres and to compute tempera- 
tures using the appropriate partition functions. It is found that the com- 
puter values closely follow the temperature profile of the chosen model. It 
is further found that the required measurable factors, namely the power ratios 
and their spatial derivatives, exhibit an appreciable sensitivity to tempera- 
ture variations, even over relatively short path increments and for relatively 
weak absorptions. 

These calculations were based on the absorption lines in the 1-0 band, 

3^- ly + transition, of the 0^ molecule. Some of these lines lie in the 
tuning range of the ruby laser. Line parameters were estimated from data ob- 
tained from R. K. Long (Ohio State University). 

Laser emissions were assumed monochromatic and simultaneously occurring 
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at three wavelengths , two of which were centered on absorption lines on oppo- 
site sides of one branch of the band and the other lying between lines. The 
choice of lines is based on the fact the lines on one side of a branch in- 
crease in strength while those on the opposite side decrease as temperature is 
varied. In this way the sensitivity of the measurement is optimized. The 
simultaneity and superposition of the pulses allow the assumption of identical 
atmosphere paths in treating their propagation characteristics, and thus elim- 
inate pulse-to-pulse variability in the returns and allows real-time analog 
processing methods to be considered. 

Another feature of the method is that the measurement of temperature is 
an absolute rather than a relative one and does not require any tie-on-mea- 
surement at some reference point along the path. 
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RAMAN SCATTERING FROM ATMOSPHERIC NITROGEN IN THE STRATOSPHERE 

*by 

M. J. Garvey & G. S. Kent 

Department of Physics, University of the West Indies 

ABSTRACT 

The Mark II laser radar system at Kingston, Jamaica, has been used to 
make observations on the Raman shifted line from atmospheric nitrogen at 
828.5 nm. The size of the system makes it possible to detect signals from 
heights of up to ho kilometres. The effects of aerosol scattering observed 
using a single wavelength are almost eliminated, and a profile of nitrogen 
density may be obtained. Assuming a constant mixing ratio, this may be in- 
terpreted as a profile of atmospheric density whose accuracy is comparable to 
that obtained from routine meteorological soundings* 

In order to obtain an accurate profile several interfering effects have 
had to be examined and, where necessary, eliminated. These include: 

1) Fluorescence in optical components 

2) Leakage of signal at 69^.3 nm. 

3) Overload effects and non-linearities in the receiving and counting 
electronics . 

Most of these effects have been carefully examined and comparisons are 
being made between the observed atmospheric density profiles and local mete- 
orological radio-sonde measurements. Good agreement has been obtained over 
the region of overlap ( 15 - 30 Km), discrepancies being of the same order as 
the experimental accuracy (1-10$), depending on height and length of period of 
observation. 
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LASER RADAR MEASUREMENTS OF. ATMOSPHERIC POTASSIUM 

by 

F. Felix, + W. Keenliside ,+ G. S. Kentt and M. C. W. Sandford* 

ABSTRACT 


A dye laser capable of transmitting in the near infra red region of the 

spectrum has been constructed to be used in conjunction with the large Mark II 

laser system at present in existence at Kingston, Jamaica, 

Preliminary measurements have been obtained of concentration of atomic 

potassium in the 70-100 km region of the atmosphere. The data indicates the 

likelihood of a double peak in the height distribution. The lower peak, which 

is the larger, is at a height of about 82 km, the upper peak is at a height 

of 9^ kin* Although an exact value for the scattering cross-section has not 

been obtained, a reasonable approximation of this parameter yields a value of 
11 -2 

about 1-15 x 10 m for the column density of atomic potassium, which is in 
agreement with other data. 


t Department of Physics , University of the West Indies 
* Radio & Space Research Station, Ditton Park, Slough, England 
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FEASIBILITY OF REMOTE MEASUREMENT OF S0 2 AND NO 
IN SMOKE STACK PLUMES 

by 

E. K. Proctor* s M. L. Wrightf* Graham Black** 


ABSTRACT 


Enforcement of pollution abatement regulations can* in many instances* 
be abetted by use of portable remote monitoring instruments that do not re- 
quire cooperation by emitters or access to their facilities. Two of the most 
important pollutant emissions from fixed sources are SO^ and NO* especially 
from coal-fired power plants. Both of these gases have absorption and 
fluorescence spectra in the eye-safe near-ultraviolet region of the spectrum 
that offer potential mechanisms for remote measurement with lidar equipment. 
Prospective ways of realizing these potentials are analyzed in this paper. 

The effectiveness of the ozone layer in reducing the daytime background noise 
level is also examined. The problem of achieving adequate penetration of 
plumes by probing signals at the concentration levels of interest is addressed* 
and a promising way of overcoming this difficult penetration problem is 
described. 


* Radio Physics Laboratory 
t Systems Techniques Laboratory 

**Molecular Physics Group 



DETECTION OF SO AND N0 o IN STACK PLUME BY RAMAN SCATTERING 
2 2 

AND FLUORESCENCE 

S. Nakahara, K. Ito and S. Ito 
Mitsubishi Electric Corporation, Kamakura Works 

325 Kamimachiya , Kamakura , Japan 
ABSTRACT 

We have been studying laser-Raman radar which can be used as the remote 

detector of S0 o concentration in the stack plume of boiler exhaust gas , and 

2 l) 

some results have been published 

In this paper , we report the interference of fluorescence against 
SO^ Raman scattering and the measuring method of SO^ and concentration. 

In a stack plume, high density dust and high concentration CO^ are in- 
cluded, therefore very strong Mie back-scattering and CO^ Raman scattering 
are observed. The separation of these scattering signals from SO^ Raman sig- 
nal was the first problem for the laser-Raman radar. But, this problem was 
solved by using the filter which have high resolving power. 

It is well known that NO^ can be excited to emit fluorescence of broad 
spectrum by blue-green light. The light source of the laser-Raman radar is 
SH of Nd:YAG laser, and boiler exhaust gas includes several tens ppm of WO^ , 
then the interference of fluorescence brings the error to SO^ measurement 

The rejection of the interference can be achieved by the subtraction of 
the fluorescence contribution from detected signal at SO^ Raman scatter- 
ing wavelength. The NO^ fluorescence contribution can be measured by two 
methods. The first is to convert the fluorescence intensity which is 
measured at a different wavelength from SO^ Raman line into the one at SO^ 
Raman line. The second is to convert the detected signal intensity, which is 
obtained when the range gate is set just behind the plume, into the one ob- 
tained at the range of the plume using the difference of the time dependence 
between SO^ Raman scattering and fluorescence. NO^ fluorescence has life 
time of about 300nS, while Raman scattering has none. 

By either of the two methods, the contribution of fluorescence is de 
termined. Then, we can measure the and SO^ concentration. 
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The minimum detectable concentration and the experimental results of the 
remote sensing of SO^ and RO^ in stack plume will be discussed in detail. 

Reference 

l) S. Rakahara et al. Opto-Electronics 4(1972) 169VL77 • 
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LIDAR OBSERVATIONS OF RAMAN SCATTERING FROM S0 2 
IN A POWER PLANT STACK PLUME 
By M. L. Brumfield, S. H. Melfi , and R. W. Storey, Jr. 

ABSTRACT 

LIDAR techniques have been successfully applied to the detection of the 
Raman backscatter from SO^ in the plume of a 200 megawatt coal -burning elec- 
trical-generating plant from a distance of 210 meters. The LIDAR system used 
consists of a 6l cm diameter, f A Newtonian telescope and 1.0 - 1.5 joules- 
per-pu'lse, 1 pulse-per-second ruby laser. Narrow band interference filters 
are used to select the T546A vibrational line of SO^, The signal from a 
photomultiplier tube was sequentially applied to each 25^ nsec wide channel of 
a 15 -channel photon counting system, resulting in a direct correlation between 
channel number and range increment. Photon counts were accumulated from the 
backscatter of a number of laser pulses (typically 50 or 100), and the accumu- 
lated counts per channel printed on paper tape. 

One sequence of measurements was made during a two-hour period while the 
plant electrical output was being reduced by approximately 50%. Although the 
Raman system had not been quantitatively calibrated, the LIDAR data correlated' 
well with the varying plant electrical output. N^ scattering observations 
were also made and an approximate quantitative SO^ concentration obtained by 
ratioing the SO^ data to N^ data. This ratio compared well to the in-situ 
measurements made during the same period by Environmental Protection Agency 
sampling instruments. 
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RAMAN LIDAR MEASUREMENT OF AIRCRAFT 
TURBINE ENGINE EXHAUST EMISSIONS 

Donald A. Leonard 
Avco Everett Research Laboratory 
Everett, Massachusetts 021^9 

ABSTRACT 

A program is currently underway to demonstrate the use of laser Raman 
backscattering as a means to measure the Composition and temperature of air- 
craft engine exhaust emissions. 

The physical property of turbine exhausts that most significantly influ- 
ences Raman spectroscopic analysis is the relatively low pollutant concentra- 
tions in combination with elevated temperatures such that significant upper 
level excitations can occur with subsequent broadening and overlapping of the 
Raman spectra. These effects must be correctly taken into account if Raman 
scattering is to be used for quantitative gas analysis at elevated temperatures. 

An experimental lidar unit was constructed using a 0.5 watt average power 
3371 & pulsed nitrogen laser, a 24-inch diameter Dall-Kirkham Cassegrainian-type 
receiver telescope and a 1 meter, double, scanning spectrometer for spectral 
analysis. The entire system operation is controlled by computer, including 
calibration, data collection and data analysis. 

Preliminary experimental results will be presented. 
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A METHOD OF MONITORING NON-RESONANT RAMAN LIDAR 
RETURNS DURING DAYLIGHT HOURS 

John Cooney 
Drexel University 

Physics and Atmospheric Science Department 
ABSTRACT 

A method for using non-resonant Raman lidar during daytime hours is pre- 
sented. 

A dual (or quadruple) arrangement of electro-optical channels is set up 
in the lidar receiver. Each channel contains its own optical spectral re- 
sponse as determined hy an appropriate set of interference filters. The in- 
coming signal is split and fed into two channels , filtered (simultaneously) 
and then inserted into the input terminals of an operational amplifier. The 
difference signal is selected and amplified whereas the common portion of the 

signal is highly attenuated (Common Mode Rejection Ratio = 75-80db). For the 

2 3 

present hardware, daylight intensities 10 -10 times Raman signals are re- 
ceived simultaneously. Mean values of intensities of daylight signals which 
give rise to a constant differences merely shift the value of the measured 
variable at the ground. These mean value signals have their common part 
highly attenuated. 

The differences of the fluctuation portions of the daylight signal repre- 
sents a noise signal. Identical (space-time) optical paths for each component 
of the daylight signal produce the same fluctuational spectra in each channel. 
Hence noise levels, much less than signal levels, are anticipated. 

First a N^-H^) signal pair is differenced. Then an N^-O^ 
enced and used for instrumental normalization purposes. 

Signal levels to be encountered are quite well known as a result of prior 
field work. 

Improved optical design of laser system can reduce mean value of daylight 

p 

intensities by 10 . Thus the mean daylight intensity can be set equal to the 

expected signal levels by more advanced optical design. 

Performance characteristics of photomultiplier pairs are discussed in 
this connection. 
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LIDAR DEVELOPMENT AT SRI - THE FIRST DECADE 


R. T. H. COLLIS 

Atmospheric Sciences Laboratory 

Stanford Research Institute, Menlo Park, California 

ABSTRACT 

On 22 July 1963 the first lidar observations of the lower atmosphere 
were made at SRI with a pulsed ruby system developed by a team led by the 
late Myron G. H. Ligda. Since that time SRI has carried out a continuous 
program of exploration and development of the technique, primarily related to 
applications in atmospheric research in the troposphere. 

In this frankly personal retrospective, some of the highlights of this 
ten years are reviewed, both in terms of progress made and difficulties ex- 
perienced. Topics discussed will include the technological aspects of the 
lidar systems used, the range of applications identified and explored and the 
various forms of information recovery and display that have been developed. 
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LIDAR SCATTERING IN THE TROPOSPHERE 


A. I. Carswell, J. D. Houston, 

W. R. McNeil, S. R. Pal 

Centre for Research in Experimental Space Science 
and Department of Physics 

York University, 4700 Keele St., Toronto, Canada 

M3J 1P3 ' 

ABSTRACT 

This paper describes measurements of the properties of the backscattered 
signals from a ruby laser radar obtained under a variety of atmospheric condi- 
tions. The lidar system constructed for these measurements employs the funda- 
mental (694.3nm) and second harmonic (34-7 . 2nm) ruby wavelengths. A Pockel’s 
cell and polarizer in the cavity are used to provide a linearly polarized 
Q-switched output with energies up to 3 Joules at the fundamental. A KDP 
doubler external to the cavity provides the second harmonic. The pulse dura- 
tion is approximately 15 ns and the pulse repetition rate is about 10 pp min. 
The transmitter collimator provides an output beam 5 cm in diameter with a 
divergence less than one milliradian. 

Pour optical receivers are equipped with rotatable polarizers and retar- 
dation plates , interchangeable narrow band filters and photomultipliers to 
detect the backscatter signal. The main receiver channel is a 20 cm diameter 
Newtonian telescope and the other three are refractive telescopes with 7*5 cm 
apertures. With this system it is possible to measure the complete polariza- 
tion properties of the backscattered signal at either wavelength with a single 
transmitted pulse. This feature permits the study of rapidly changing atmos- 
pheric properties whose polarization properties could not be obtained by se- 
quential firings. 

A parameter of particular interest is the depolarization, 6, which is de- 
termined from the relation 6 = P /P„ where P and P n denote the backscattered 

r Z r x, 

power polarized perpendicular and parallel to the trnasmitter polarization. 
This parameter is very easily measured and is found to vary greatly with at- 
mospheric conditions. Measurements with this system have been in progress for 
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about two years and provide considerable information on the variation of the 
backscattered signal polarization for a variety of meteorological conditions 
at altitudes up to about 10km. For "clear” (urban) air 8 is typically of the 
order of 0.2 but it can vary between 0.05 and 0.9 while well developed hazes 
have 6 generally smaller than 0.1. For fog, 8 values at both wavelengths are 
around 0.1 and higher. Snow gives values of 8 . of the order of unity. Clouds 
exhibit depolarizations covering the whole range from about zero to unity and 
in many cases the presence of multiple scattering is clearly seen in the in- 
crease of depolarization with penetration depth. Cloud measurements show con- 
siderable variation of 6 with altitude and the results indicate that polariza- 
tion signatures could be useful for cloud characterization and classification. 
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LIDAR MEASUREMENTS OF THE THERMIC STRUCTURE IN THE TROPOSPHERE 

Ch. Werner 

German Establishment for Aeronautical and Space Research (DFVLR) 
Institute for Atmospheric Physics 

D-8031 Oberpf af f enhof en , Germany 
ABSTRACT 

Lidar measurements by UTHE (1972) and similar measurements with acoustic 
radar by HALL (1972) have shown, that short time variations in the troposphere 
can be determined. 

Aerosol backscatter measurements give the vertical structure. A series 
of such results shows the time variation of the tropospheric structure. The 
Figure is an example for a measuring period during half an hour in summer 1972 
(l pulse per minute). These variations were produced by humidity variations. 
For example, at humidity variations from 75$ to 77$ aerosol particles with 1 
micron radius change their radii during 50 to 100 milliseconds. By lidar 
technique it is possible to determine this effect qualitatively. With a pulse 
repetition rate of 1 per second we plan to investigate the thermic structure 
at different weather conditions. 

By simultaneous measurements of the nitrogen Raman component it would be 
possible to compare the short time variations by temperature and humidity. 

The influence of wind can be determined by scanning a defined region (region 
or space). 

The meteorologist is interested in results of short time variations in 
the troposphere (microscale). By this remote sensing method he can get 
these values continuously. 


UTHE, E. E. 1972: 1+th laser radar conference, Tucson 
HALL, F. F. 1972: Short Course in Remote Sensing, Boulder 
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Time variation of the tropospheric structure 


LIDAR MEASUREMENTS OF ON-SHORE WIND DIFFUSION 


R. M. Brown, P. Michael and G. Raynor 
Brookhaven National Laboratory 
ABSTRACT 

The concept to place electric power generating stations on the oceans 
off the coast of the United States has instilled new efforts in research for 
improved understanding of the diffusion properties of the atmosphere in the 
ocean-air interface. 

The Atomic Energy Commission has instigated a program by the Meteorology 
Group at Brookhaven National Laboratory to investigate the low level, on-shore 
wind systems that dominate many of the coastal regions. Analytical techniques 
and specialized instrumentation from previous studies at Brookhaven are being 
used in this new program. The Brookhaven Lidar system is used to measure some 
of the physical properties of the oil-fog plume originating from a portable 
smoke generator on a boat off the coast. The oil-fog plume is used as a 
tracer which can be observed, photographed and measured to determine the dif- 
fusive power of the atmosphere associated with the ocean-air interface and 
the discontinuities found in the ocean-land boundary. 

This paper will describe the program rather briefly and the oil-fog 
scattering measurements that have been made with the Lidar system. 


This abstract submitted for presentation at the Fifth Conference on Lasar 
Radar Studies of the Atmosphere, June k-6 , 1973, Williamsburg, Va. 
Sponsored by 1. Group on Laser Atmospheric Probing 

2. American Meteorology Society 

3. Optical Society of America 
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ESTIMATING VERTICAL DIFFUSION COEFFICIENTS BY LIDAR 

Walter M. Culkowski 
Searle D. Swisher 

Atmospheric Turbulence & Diffusion Laboratory 

National Oceanic & Atmospheric Administration 

Post Office Box E 

Oak Ridge, Tennessee 37830 
ABSTRACT 

The Atmospheric Turbulence and Diffusion Laboratory at Oak Ridge, 

Tennessee has been conducting routine probing of the lower troposphere and 

comparing the results with those obtained with turbidity photometers and a 

distant suspended particulate station* The change in scale height, K^/v , 

with time permits the vertical turbulence coefficient to be estimated if 

V is known or assumed. Extremely high monthly correlations of turbidity vs. 
s 

the log of backscatte.r at 100 meters have been obtained. In addition, high 
correlations of suspended particulate matter at Chattanooga and Oak Ridge sug- 
gest that the bulk of particulate matter is of natural, rather than industrial, 
origin. 
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ATMOSPHERIC AEROSOL AND THERMAL STRUCTURE 
IN THE BOUNDARY LAYER OVER THE LOS ANGELES BASIN* 


Warren. B. Johnson, Ph.D 
Atmospheric Sciences Laboratory 
Stanford Research Institute 
Menlo Park, California 9^-025 
ABSTRACT 

A field study using a mobile lidar was recently conducted in the 
L. A. Basin, California, to (l) examine the relationship between the 
vertical aerosol and the thermal structure, and (2) map the vertical 
aerosol structure in the atmospheric boundary layer over the basin. 

These data are needed for use in the development of a mixing-depth 
submodel required for photochemical air quality simulation models. 

Toward these ends, a series of lidar aerosol measurements in conjunction 
with balloon and aircraft temperature soundings were taken at a site in 
El Monte, and in a mobile mode along a 90-mile freeway loop between 
El Monte, Santa Monica, and Long Beach. The lidar data are presented in 
the form of time-height and distance-height cross sections. The results 
indicate that, although aerosol concentrations are frequently present 
above the base of the marine inversion, these are generally in stratified 
layers in contrast to the more uniform nature of the lower convective layer, 
permitting the mixing depth to be distinguished on this basis. The lidar- 
derived mixing depths are well correlated (within 100 m) with daytime 
temperature inversions. Other significant features shown by the lidar 
data include large Basin-wide mixing-depth variations, waves with 
amplitudes of 200-300 m and wavelengths of 1000-1500 m on the lower aerosol 
layer, and apparent aerosol "chimneys" with overrunning in the vicinity of 
convergence zones . 


This work was performed while the author was affiliated with the Meteorology 
Laboratory, U.S. Environmental Protection Agency, Research Triangle Park, N.C. 
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LASER SOUNDING OF INDUSTRIAL MIST 


V. E. Zuev 9 B. V. Kaul, I. V. Samokhalov 

Institute of Atmospheric Optics of the Siberian 
Branch of the U.S.S.R. Academy of Sciences, 

Tomsk, U.S.S.R. 

L. S. Ivlev, K. Ya. Kondratiev 

Leningrad State University, Leningrad, U.S.S.R. 

ABSTRACT 

A ruby lidar has been installed at about 2 km from the pollution sources. 

The profiles of the backscattering factor c and attenuation factor 
a are obtained by the probing in different directions. The absolute laser 
calibration has been realized that permits to solve a laser radar equation 
using layered restitution of profiles a and a. The value of a phase 
function y 1 in l 80 ° direction is estimated in the surface boundary layer 

assuming the atmosphere being homogeneous in a horizontal plane. In this case 
a logarithm of the signal value multiplied by distance square is a straight 
line according to the slope of which a can be found and the laser radar 
equation can be solved. 

Thus estimated values of y oscillate about average value of 0.033 

with the root-mean-square deviation +_ 0.0105, An average value of a phase 

function is used for solving the laser radar equation when it is impossible to 

estimate y . 

'tt 

The solutions obtained allow to derive the field of 0 values and coin- 
'll 

pare these data with dustiness estimated by direct sampling. A comparison 
proves to be possible with an error up to 70 per cent. Feasible reasons for 
obtaining poor agreement of the optical measurements with the results of 
sampling are discussed. 

Some attempts have been taken up to restore the profile at a along 
vertical and inclined routes. An essential inhomogeneity of the atmosphere 
is observed in the horizontal direction up to 1800 m height. This makes 
determination of vertical transparency by probing along different zenith 
directions to be of little use under the industrial city conditions. Estimate 
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of the concentration gradient of contaminating impurity proves to be possible 
at different altitudes. 

The layers of increased turbidity are observed at 1200-1700 m height. 

The statistical relation between the existence of such layers and the dustiness 
of the surface boundary layer is found. The coefficient of correlation between 
the value of at 1500 m level and the value of dustiness in the boundary 

layer is estimated to be 0 . 72 . 
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TROPOSPHERIC TRMSMISSIVITY MEASUREMENTS USING THE 
RAMAN NITROGEN LIDAR TECHNIQUE 

toy 

M. P. McCormick, G. B. Northam, and W. H. Fuller 
NASA - Langley Research Center 
ABSTRACT 

LIDAR measurements in Azusa, California, during October 1972, were made 
in which the backscattered Raman-shifted nitrogen return was ratioed at dif- 
ferent altitudes in order to obtain transmissivity. Rawinsonde data from 
nearby El Monte were used to determine the temperature and nitrogen number 
density altitude profiles. 

These data and other meteorological data are compared to the vertical 
aerosol and transmissivity structure determined by LIDAR. Also data analysis 
techniques are shown for obtaining q (transmissivity) and 0 (attenuation 
coefficient) as a function of altitude. 
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A COMPARISON OF ATMOSPHERIC STRUCTURE AS 


OBSERVED WITH LIDAR AND ACOUSTIC SOUNDER TECHNIQUES 

Edward E . Uthe 
Stanford Research Institute 
Menlo Park, California 

Neil A. Shaw 

Argonne National Laboratories 

Argonne, Illinois 
ABSTRACT 

During the METROMEX* 1972 field program, Stanford Research Institute (SRI) 
and Argonne National Laboratories (ANL) conducted a joint lidar and acoustic 
sounder data collection program in a manner that readily allows comparison of 
atmospheric structure observed with these instruments. A preliminary analysis 
of the data has been completed and it is the purpose of this paper to present 
these results. 

The SRI Mark IX lidar uses a Pockels cell Q-switched ruby laser capable 
at firing at 60 pulses/min with an energy output of approximately 1.0 joules/ 
pulse. A 2-inch transmitter lens is coaxial with a 6-inch Newtonian receiver. 
Lidar backscatter signatures were electronically corrected for the inverse 
range-squared dependence, logarithmically amplified, and stored on a magnetic 
video-disc recorder. Each of the 1000 tracks of the disc can contain up to l6o 
lidar signatures and these signals can be fed through electronics that produce 
intensity-modulated cross-sections of aerosol structure. 

The ANL acoustic sounder emits sound pulses of approximately 100 ms in 
duration, 100 watts (RMS) in power, and near 1000 Hz in frequency. The 
transmitter and receiver uses a five-foot diameter parabolic dish with the 
transmitter/receiver transducer mounted at the focus. During the field 

*METR0p olitan Meteorological Experiment — an investigation of weather modifica- 
tion by urban effects. The SRI program was supported by NSF Grant GI-3^770 
(Environmental Systems and Resources, Weather Modification Program). 
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program, a sound pulse was emitted every 12 s and the detected backscattered 
energy from each pulse was filtered, logarithmically amplified, corrected for 
the inverse range- squared dependence, and displayed on a facsimile- type 
recorder. Hence, the lidar and sounder returns were similarly processed and 
displayed and, therefore, can be readily compared. 

The displayed returns (S) can be related to the single-scattering equation: 

S = log = log (k 2 6T 2 ) 

Where P = power returned 

R - range 

K^, = instrumentation constants 

E = transmitted energy 

3 - volume backscatter coefficient 

2 

T = round-trip transmission. 

2 

For the lidar, 3 and T are principally determined by the optical properties 
of the scattering particles and the particle concentrations. For the sounder, 

o 

3 is determined by temperature irregularities (fluctuations) and T is 
determined by scattering from temperature and wind fluctuations and absorption 
by water vapor. Hence, the backscatter returns are related to entirely different 
atmospheric parameters; however, they may be spatially and temporally corre- 
lated because of similar dependence on various atmospheric dynamic and physical 
processes . 

Figure 1 presents an example of the data collected with the lidar and the 
sounder separated by approximately 100-yards and both pointing vertically 
upward with a lidar beamwidth of one milliradian and a sounder beamwidth of 
120. The lidar firing rate was 3 pulses/min. Increased signal return in these 
presentations is represented as increased brightness for the lidar data and 
decreased brightness for the sounder data. (Note that there is a factor of two 
difference in the vertical scales of the two records.) The dark streaks in the 
sounder record are a result of background noise sources, primarily jet aircraft. 

The early morning mixing depth and its gradual increase by convection as 
determined from both data types are well correlated. However, the sounder 
record does not give this information after 1230 CDT. This is a result of low 
signal-to-noise ratios at heights greater than 3000-feet and because vertical 
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air motions associated with active convection tend to give returns containing 
little vertical structure. The turbulent cells are easily identified on the 
lidar presentation (1330-1500 CDT). During this time period, near-surface wind 
sensors recorded light and variable winds. The wind speed increased and shifted 
in direction at 1530-1620 CDT. These events are recorded on the lidar data as 
sudden decreases in the near surface aerosol concentrations and on the sounder 
data as sharp vertical discontinuities. 

A cellular-structured air parcel passed over the area between 1800 and 
1900 CDT that lifted the atmospheric layers. This perturbation introduced 
short-period atmospheric waves that were detected by both indirect sensors. 

The sounder record during this time period indicates that a new low-level 
mixing layer was being established. The sounder normally detects the 
reestablishment of the mixing layer before the lidar, since particulate gradients 
normally result from a trapping of particles by thermally stable stratifications. 

Other data to be presented illustrates complex day-to-day variations of 
atmospheric structure and provides information on how lidar and sounders may 
effectively be used in atmospheric research and operational programs. 
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OBSERVATIONS OF INVISIBLE PLUME BY SHGed YAG LASER RADAR 


H. Tanizaki, F. Ito , T. Yakuo, I. Yamashita, T. Hosokawa 
Nippon Electric Company 5 Fuchu Plant 
1-10 5 Nisshin Cho , Fuchu City, 183 Japan 
(0423)64-1111 

T . Ikeno 

Japan Weather Association 
ABSTRACT 

The enormous* efforts ^to reduce exhausted materials contained in a plume 
which cause air pollution has been taken into plant operations. As a result , 
exhausted plume became invisible to naked eyes and new observation facilities 
are needed to monitor the diffusion of plume. The following is the report on 
the experimental observation to confirm the suitability of SHGed YAG Laser 
Radar for invisible plume monitoring. 

A certain Kombinat area near Tokyo was chosen as a test range. The ex- 
hausted plume from ultra high stacks which we chose as a target were invisible 
to naked eyes even in clear daytime. Eight observations , three in daytime 

and five in nighttime, were performed at two locations within the Kombinat 

l) 

yard using mobile SHGed YAG Laser Radar which has 1 MW output power at 
0.53 micron, 50 pps repetition frequency and facilities of RHI 5 PPI and A scope 
display. The RHI and PPI display were recorded on 35 mm photographic films 
for further analysis. 

The photographic data of three observations, one in daytime and two in 
nighttime, were chosen and analysed. 

Fig. 1 and fig. 2 show the analysed result of observation performed on 
1972, 9 >25 at 2100 - 2120 JST. Fig. 1 is plan position indication which is the 
projection of diffused volume to the ground and fig. 2 is range height indica- 

l) H. Tanizaki et al; Conference abstracts of the fourth conference on laser 
radar studies of the atmosphere (Jan. 26 - Jan. 28, 1972, Tucson, Arizona 
U.S.A. ) pp. 82 


88 



tion which is the projection of diffused volume to the vertical plane parallel 
to the wind direction. The wind direction, wind speed and temperature at the 
observation were ENE, 1 m/sec and 21 degree Centigrade respectively. 

It is obvious from fig. 1 and fig. 2 that the diffusing space expanded to 
the wind direction and the effective height of three ultra high stacks reached 
respectively 1.7 » 1.8 and 1.5 times of their physical heights at that time 
presumably due to slow wind speed. 

Thus, the suitability of SHGed YAG Laser Radar for invisible plume dif- 
fusion monitoring is confirmed by this observation as we have been technically 
anticipated. But, real time Laser Radar echo processing system must be 
developed before practical use of Laser Radar as a monitoring system. 
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DETECTION AND ESTIMATION OF ATMOSPHERIC 


TURBIDITY FROM THE POLARIZATION OF SCATTERED RADIATION 

C. R. Nagaraj a Rao 

Atmospheric Optics and Radiation Lahoratory 
Department of Meteorology 
University of California 
Los Angeles 

ABSTRACT 

Estimates of the 'effective* turbidity of the atmosphere due to particu- 
late pollution have been made from an interpretation of the results of ground- 
based and high altitude balloon measurements of the polarization of scattered 
radiation, obtained over the past five years, in the light of numerical compu- 
tations of the polarization of emergent radiation in models of turbid atmos- 
pheres. Ground-based measurements were made at Los Angeles where the occur- 
rence of low level pollution is frequent and well defined. The balloon mea- 
surements were made from altitudes of 28-30 km over southwestern New Mexico. 
Programmed, servo-governed photoelectric polarimeters were used in these 
measurements which were confined to the near ultraviolet and visible regions 
of the spectrum. 

The polarization of emergent radiation was computed in two simple, singu- 
lar but possible models of a turbid atmosphere based on the solution of the 
radiative transfer problem. The variable model parameters are the amount and 
type of aerosols and the reflective properties of the ground. Computations 
have been made for two different types of polydisperse aerosols with sizes 
governed by the modified gamma and power law distributions and with different 
refractive indices. The ground has been assumed to be a Lambertian reflector. 

The optimal values of the aerosol optical thickness which yield good cor- 
respondence between theory and experimen 4- are found to depend upon the phase 
function asymmetry factor associated with the type of aerosols used In the 
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model computations. The turbidity factor, defined as the ratio of the total 
optical thickness of the atmosphere to the molecular optical thickness, varies 
from 2 to 5*5 over the spectral region of interest on normal turbid days. The 
similar effects that the aerosols and the Lambertian ground have on the emer- 
gent radiation introduces a certain amount of ambiguity in the interpretation 
of the experimental results. This ambiguity will be discussed in terms of a 
model parameter sensitivity analysis. 
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ON THE INVERSION OF LIGHT SCATTERING DATA 


INTO INFORMATION ON AEROSOL PROPERTIES 

by 

A. Cohen E. Westwater, R. E. Cupp and V* E. Derr 

Wave Propagation Lab., NOAA-ERL, Boulder , Colo. 

ABSTRACT 

The Backus-Gilbert inversion technique was applied to scattering data 

and calculated examples were given to ascertain its ability to recover the 

1 

aerosol properties with no a-priori assumptions. 

In order to cheque this inversion procedure as well as other techniques 

(such as Daves' and Herman's et al. ) , a controlled experiment was performed. 

The experiment dealt with spherical particles the Mie scattering curves 

of which was measured by use of a dye-laser. These curves were used to ac- 

2 

curately determine the refractive index of the particles. 

Mixtures of three different known sizes were then prepared and the 
scattering intensity vs. wavelength was measured at 90° scattering angle. 

The mixtures contained also various known relative concentrations of the 
three sizes. 

The analysis of the measured results will be shown and discussed. 


1. E. Westwater and A. Cohen Appl. Opt. June 1973* 

2. A. Cohen, V. E. Derr, R. E. Cupp and T. McNice Appl. Opt. May 1973. 

+ Now with the Dept, of Atmos. Sciences, the Hebrew University of 
J erus alem , Israel . 
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THE AEROSPACE CORPORATION AEROSOL SCATTERING STUDIES 

Richard T. Hall, R. Douglas Rawcliffe, 
and David C. Cart-wright 
ABSTRACT 

An experimental apparatus for measuring the differential scattering cross 
section of micron-sized aerosol particles has been designed and constructed 
at the Aerospace Corporation. The apparatus uses the 10.6 ym laser transi- 
tion of CO^ for illumination of the aerosol particles. Two mercury-doped 
germanium, liquid helium cooled detectors are used to measure the scattered 
10.6 ym radiation. One detector is mounted at right angles to the laser 
beam while the other is moveable between the forward scattering and 
backwards scattering positions. 

The apparatus has been used to measure the differential scattering 
cross sections of 9 to 20 y diameter carbon particles and 2 to 15y diameter 
aluminum particles. Calculations of the differential scattering cross 
section as predicted from Mie theory have been carried out to provide a 
comparison between theory and experiment. 
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SIMULTANEOUS RED - BLUE LIDAR 
AND AIRBORNE IMP ACTOR MEASUREMENTS 


•by 

M. P. McCormick 9 I* H.. Blifford° 

* o 

W. H.. Fuller 5 and G. W. Grams 

ABSTRACT 

Simultaneous two-color (0.6943 micro meters and 0.3472 micro meters) 
LIDAR measurements were made in the troposphere and lower stratosphere over 
Boulder 5 Colorado during March 1973* In addition, on the evening of March 26 , 
airborne single-stage impactor measurement; were made at four altitudes — 
10,500, 25 ,000, 33,000 and 43,000 feet MSL. These data were integrated at 
constant altitude for 15, 45, 45, and 60 minutes respectively. 

The LIDAR data were taken with Langley *s 48" LIDAR using a dichroic 
beamsplitter to separate the return at 0*6943 micro meters and 0.3472 micro 
meters. The analog waveforms for both colors were digitized simultaneously; 
one on an RCAR data acquisition system and the other on the 48" Langley data 
acquisition system. A discussion of the preliminary results from these 
measurements will be presented. 


*NASA -• Langley Research. Center 
National Center for Atmospheric Research. 
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MULTIPLE SCATTERING MEASUREMENTS AS A 


FUNCTION OF WAVELENGTH BY USE OF A DYE-LASER 

Dy 

Ariel Cohen,+ Richard E. Cupp and Vernon E. Derr 

Wave Propagation Laboratory 

NOAA-ERL, Boulder, Colorado 
ABSTRACT 

Dow latex monosized spheres were used as scatteres in a laboratory sys- 
tem containing a dye-laser as the light source (the system has been described 
1 

elsewhere ). Eight samples of different particle concentrations (1:2:4:8:16: 
32:64:128) were prepared and the optical depth of each was measured. The two 
lowest concentrations gave the same relative Mie Scattering intensity curve as 
a function of wavelength as is expected when merely single scattering occurs, 
(the absolute intensity was doubled at all wavelengths). The shape of the 

\ 

curve started to change with the further increase in the number density of the 
particles. The main characteristics of the changes was the gradually disap- 
pearance of the minima and for higher concentrations - also the maxima. 

Special measurements were made to ensure that the beam path and the sample 
volume were big enough to include most possible multiple scattering effects 
influencing the scattering intensity. Thus, this experiment provides (within 
a specified degree of accuracy) the multiple scattering Mie curves vs. optical 
depth for a range of ~ 1000 5L ■>£; 


1. A„ Cohen, V. E. Derr- T. McNice and R. E. Cupp, Mie Scattering.,... 
Appl. Opt. May, 1973. 

+ Now with the Dept, of Atmospheric Sciences, The Hebrew University of 
Jerusalem, Israel. 
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AIRBORNE LASER POLAR NEPHELOMETER 


*y 


Gerald W. Grams 

National Center for Atmospheric Research"^" 

Boulder , Colorado 
ABSTRACT 

A polar nephelometer has been developed at NCAR to measure the angular 
variation of the intensity of light scattered by air molecules and particles* 
The system has been designed for airborne measurements using outside air 
ducted through a 5-cm diameter airflow tube; the sample volume is that which 
is common to the intersection of a collimated source beam and the detector 
field of view within the airflow tube* The source is a linearly polarized 
helium-neon laser beam* The optical system defines a collimated field-of- 
view (0*5° half-angle) through a series of diaphragms located behind a 172-mm 
focal length objective lens. A photomultiplier tube is located immediately 
behind an aperture in the focal plane of the objective lens. The laser beam 
is mechanically chopped (on-off) at a rate of 5 Hz ; a two-channel pulse 
counter 9 synchronized to the laser output, measures the photomultiplier pulse 
rate with the light beam both on and off. The difference in these measured 
pulse rates is directly proportional to the intensity of the scattered light 
from the volume common to the intersection of the laser beam and the detector 
field-of-view. 

Measurements can be made at scattering angles from 15° to 165° with 
reference to the direction of propagation of the light beam. Intermediate 
angles are obtained by selecting the angular increments desired between these 
extreme angles (any multiple of 0.1° can be selected for the angular increment 
5° is used in normal operation). Pulses provided by digital circuits control 
a stepping motor which sequentially rotates the detector by pre-selected 
angular increments. The synchronous photon-counting system automatically be- 
gins measurement of the scattered-light intensity immediately after the rota- 
tion to a new angle has been completed. 
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The instrument has been flown on the NASA Convair 990 airborne laboratory 
to obtain data on the complex index of refraction of atmospheric aerosols. A 
particle impaction device is operated simultaneously to collect particles from 
the same airflow tube used to make the scattered-light measurements. A size 
distribution function is obtained by analysis of the particles collected by 
the impaction device. Calculated values of the angular variation of the 
scattered-light intensity are obtained by applying Mie scattering theory to 
the observed size distribution function and assuming different values of the 
complex index of refraction of the particles. The calculated values are then 
compared with data on the actual variation of the scattered-light intensity 
obtained with the polar nephelometer. The most probable value of the complex 
refractive index is that which provides the best fit between the experimental 
light scattering data and the parameters calculated from the observed size 
distribution function. 


1 The National Center for Atmospheric Research is sponsored by the National 
Science Foundation. 
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OPTICAL MEASUREMENT OF ATMOSPHERIC AEROSOL 


“by 

A. P. Waggoner 
and 

R* J. Char Ison 
ABSTRACT 

Nephelometric measurements of extinction due to aerosol scatter are found 
to he consistent with solar and stellar extinction measurements hy Voltz photo- 
meter and astronomical telescope respectively fall measurements at or corrected 
to 500 nm] . Sites without obvious contamination by man or, regional vegetation 
such as Mauna Loa and Mauna Kea, Hawaii have optical aerosol scatter and ex- 
tinction coefficients that average 0.2 or less of Rayleigh scatter and extinc- 
tion. Average stellar aerosol extinction values of less than 0.31 Rayleigh 
extinction have been measured at Cerro Tololo, Chile; Cape of Good Hope, South 
Africa; Kitt Peak, Arizona. Solar extinction measurements at some hO stations 
at rural and -urban locations in the U.S. show that at every location there are 
clean periods when the measured aerosol extinction is less than 0.5 Rayleigh, 
and the average aerosol extinction is only 1.2 times Rayleigh 4 over the central 
western states during continental polar synoptic conditions. 

This collected data, from measurements by our group and others, is not 
consistent with either extinction calculated from searchlight measurements 
or the model atmospheres of Elterman. 


Both authors are on the staff of the University of Washington. 
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NUMERICAL EXPERIMENT ON MULTIFREQUENCY LASER SOUNDING OF CLOUDS AND ESTI- 
MATION OF PARTICLE SIZE DISTRIBUTION 

V. E. Zuev, G. M. Krekov, M. M. Krekova, E. V. Makienko, 

I. E. Naats 

Institute of Atmospheric Optics of the Siberian 
Branch of the U.S.S.R. Academy of Sciences, 

Tomsk, U.S.S.R. 

Determination of particle size distribution parameters of polydisperse 
medium on the basis of optical measurements is connected with solution of 
inverse problems of scattering theory. Solution of these problems in general 
form encounters definite mathematical difficulties. However, in the case of 
cloudy media a priori information on the refractive index and the particle 
size distribution form allows practically realizable procedures of the opti- 
cal measurements inversion to be obtained. 

An effective numerical algorithm based on the method of optimal para- 
metrization is suggested for this purpose in the present paper. As a priori 
model of distribution function the modified gamma distribution is used [l], 
the parameters of which are determined from the optical measurements or 
equivalent numerical estimates, performed for a number of wavelengths. 

Necessary information on characteristics of scattered radiation is derived 
from solution of a known laser radar equation valid for the single scattering 
approximation. It is supposed that small viewing field of optical detectors 
gives sufficient guarantee of applicability of single scattering equations to 
calculate the lidar systems. However, as the estimates show, at sounding of 
optically dense atmospheric aerosols (clouds, fogs, ground mist) this question 
requires careful research. 

Accordingly, numerous calculations of space-time characteristics of back 
scattering at cloudy layer sounding were carried out. The numerical solution 
of the nonstationary transfer equation was produced by the Monte-Carlo 
technique for boundary conditions corresponding to the real schemes of sounding 
at \ = 0 . 69 ^ 3 *, 1 . 06 ; 2 . 36 ; 3 * 51 ; 5-3 and 10.6 y. . The results obtained give a 
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possibility to analyze the influence of multiple scattering depending on the 
choice of a cloud optical model, viewing field of a detector, cloud altitude, 
etc. In particular, substitution of the value of signal being detected into 
the laser radar equation (at preliminary calibration of "numerical lidar") 
allows the high-altitude profile of volume backscattering coefficient (^(H) 

"to be measured" taking into account a distorting correction of multiple 
scattering. 

Fig. 1 shows an example of such profile for a cloud model, having drop 

-3 

concentration invariable in the altitude H, N = 100 cm . The results are 
given for a scheme of monostatic sounding at A = 2.36 um; the height of the 
cloud lower edge is 1000 m. Curve 0 corresponds to the single scattering 
component, curves l-*i are obtained by treatment of total signal arriving 
to the detector, the angular aperture being equal to 3’j 20’, 1° and 5°» 
respectively. 

The values of ft^AjH) obtained as a result of numerical experiment are 
inverse! to the particle size spectrum by the method of optimal parametriza- 
tion. An algorithm of optimal parametrization is based on the approximation 
of particle size distribution, found by a certain model distribution, the 
parameters of which are determined in scanning process of function minimum, 
characterizing a degree of proximity of calculated optical characteristics and 
measured in some points of spectral interval. The parameters evaluated are: 
modal radius and distribution density. 

Special attention in this paper is given to the questions of optimal 
choice of wavelengths of optical sounding. Possibilities of numerical 
simulation of the experiment according to the known initial data ensure 
practically absolute control of efficiency of treatment procedure. 

Curve 0 in Fig. 2 illustrates coincidence of reducible density of particle 
size distribution with initial function of gamma distribution for a cloud 
model Ci[l] with an accuracy being not less than 1 per cent. In this case the 
absolute values of 3 7r (A) for the points A = 2.36; 3.51; and 5-3 Pm, being 
not distorted by multiple scattering, are used for the treatment. The 
distributions 1-5 are obtained according to the values of a total signal of 
back scattering, arriving in the aperture ^ = 1° from the altitudes 30, 60 , 

90, 120 and 150 m, respectively, measured from the cloud lower edge. Their 
deviation from the real one characterizes increased influence of multiple 
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scattering in the course of the light beam penetration into the cloud. 

The points of spectral interval X = 0.69^3, 1.06 and 10.6 ym appear to 
be less informative for the reasons of remote determination of size distribu- 
tion of cloud drops. 

REFERENCES 


1. D. Diermendjian, Electromagnetic Scattering on Spherical Polydispersions, 
Amer. Elsevier Publ. Co. Inc., N.Y., 1969* 
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ON LIDAR SOUNDING OF THE ATMOSPHERE TO ESTIMATE 
STATIC AND DYNAMIC CHARACTERISTICS OF AEROSOL INHOMOGENEITIES 

V. E. Zuev 9 G. 0. Zadde 

Institute of Atmospheric Optics of the Siberian 
Branch of the U.S.S.R. Academy of Sciences * 

Tomsk* U.S.S.R. 

V. P. Tarasenko* N. I. Yurga 
Siberian Physico-Technical Institute* Tomsk* U.S.S.R. 

ABSTRACT 

1. A possible application of intensity fluctuations of a pulse light 
signal reflected by atmospheric aerosols is analyzed by the correlation method 
to evaluate static (medium sizes* shape) and dynamic (speed and direction of 
movement* lifetime) characteristics of aerosol inhomogeneities. The aerosol 
inhomogeneities are assumed to be expanded* pressed, disintegrated and orig- 
inated constantly in accordance with random laws* the set of inhomogeneities 
as a whole traveling together with air masses and having predominant movement 
in wind direction. It is shown that the characteristics of aerosol inhomo- 
geneities considered can be expressed by the coefficients of the correlation 
function expansion 'of the reflected signal fluctuation intensity in Tailor 
series. 

2/ Correlation systems for evaluating static and dynamic characteristics 
of driving objects can be divided into two types according to the kind and 
quantity of used information; the systems with coordinates of the information 
removal "points” to be fixed in space* and the systems with a parallel simul- 
taneous information removal at discrete moments of time. 

The systems for determination of wind direction considered in [l*2] are 
the examples of the first type system. However* the operating information 
removal for two points is insufficient to estimate completely static and dynamic 
characteristics of inhomogeneities, their quantity ought to be increased up to 
three of them for two-dimensional problem and up to s four of them for three- 
dimensional problem as it is usually done in the ionospheric studies. 

The second type systems are used for the investigation of a medium shape 
and speed of the clouds according to photographs made from satellites. These 
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systems are also used for solution of navigation problems [ 3 ] • 

The use of optical quantum generators with a scanning beam is seen to 
increase greatly the working information removal in comparison with the first 
type systems. Nevertheless, scanning rate is not sufficient sometimes in order 
to consider a general picture of aerosol inhomogeneities to be stationary. In 
this connection the use of the systems of second type treatment becomes a mat- 
ter of essential difficulty. 

3* Aerosol inhomogeneities simulation has been carried out on the basis 
of the digital computer experiments with the aim of estimating static and 
dynamic characteristics of inhomogeneities by an optical beam in the atmosphere 
at different scanning procedures. 

The dependence of determination accuracy of these characteristics on the 
type of chosen laws of aerosol particle distributions in the atmosphere, the 
parameters of inhomogeneities geometry, their speed and the law of scanning 
have been obtained. 


REFERENCES 

1. V. E. Derr, C. G. Little, "A Comparison of Remote Sensing of the Clear 

Atmosphere by Optical, Radio, and Acoustic Radar Techniques", App. Opt., 
1970 , Vol. 9, No. 10. 

2. A. Cohen, V. E. Derr, R. E. Cupp, J. J. Barth, "On the Correlation of 
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A COMPARISON OF NEAR SIMULTANEOUS LIDAR RETURNS 
AND PARTICULATE COLLECTIONS ON FILTERS 
FLOWN AT SIX STRATOSPHERIC ALTITUDES 


Burton G-. Schuster, Allan L. Lazrus , 
and Frederick G. Fernald 
National Center for Atmospheric Research* 

Boulder, Colorado 

William Sedlacek and Paul Guthals 
Los Alamos Scientific Laboratory** 

Los Alamos, New Mexico 

ABSTRACT 

Collections of particulates on both Los Alamos Scientific Laboratory 
(LASL) and National Center for Atmospheric Research (NCAR) filter systems were 
made from an RB 57F aircraft flown at one tropospheric and six stratospheric 
altitudes over the Boulder, Colorado, area. This daytime flight was spanned 
by lidar returns on evenings before and after the flight. 

Scanning electron microscope examination of the LASL filters showed no 
evidence of solid particulates greater than .2 ym (the instrumental resolving 
power). Quantitative analysis of the NCAR filters yielded chemical composition 
and mass. The mass values were normalized to the total air flow through the 
filters to yield mass mixing ratios at the various altitudes. The lidar returns, 
normalized to molecular densities obtained from sonde data, were put in the form 
of particulate scattering divided by molecular scattering, i.e., an optical 
mixing ratio. A plot of the optical mixing ratio versus mass mixing ratio, in 
the stratosphere, yielded linear relationship, for five of the six data points, 
going through the origin. 

*The National Center for Atmospheric Research is sponsored by the National 
Science Foundation. 

**The Los Alamos Scientific Laboratory is operated under AEC contract by the 
University of California. 
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USE OF POLARIZATION LIDAR FOR INVESTIGATION 
OF METEOROLOGICAL FORMATIONS 

Yu.S.Balin, G.O.Zadde, V.E.Zuev, G.M.Krekov, 

G.G. Matvienko 9 A. I. Popkov, I .V.Samokhvalov, 

G. A. Titov, V . S . Shanianaev 
Institute of Atmospheric Optics of the Siberian 
Branch of the U.S.S.R, Academy of Sciences, 

Tomsk, U.S.S.R. 

ABSTRACT 

This paper presents the results of theoretical and experimental, investi- 
gations of depolarization characteristics of different meteorological 
formations. Experimental investigations are carried out with a monostatic 
lidar. The ruby laser radiation is polarized in a vertical plane. The 
radiation reflected is accepted by a lens system of 150 mm in diameter and 
a viewing angle of 12 1 and further it is divided by Wollaston prism into 
the components polarized orthogonally. In this case the principal plane of 
the prism is exposed parallel with the laser polarization plane. 

Investigations show the degree of radiation polarization, reflected 
from water clouds, to be changed within 1^0.7 (seldom up to 0.6) depending 
on their density. In most cases a signal reflected from the cloud leading 
edge is polarized completely. The time shift is observed between polarized 
and crosspolarized components of a signal, reflected from a cloud, depending 
on the density of a meteorological object. While penetrating into the cloud 
depth a degree of polarization decreases up to 0.8-0. 7* and the character of 
this decrease is different for various types of clouds. 

For crystal clouds the shift between the components of the reflected 
signal is not observed and the magnitude of polarization degree .amounts to 
O.l-rO.3 in comparison with water clouds. The polarization degree of radiation 
reflected by fog is not less than 0.6, and that in the rains of average 
intensity (about 5im/h) I s always about 1. 

The authors have suggested an algorithm of numerical solution of 
nonstationary transfer equation in the vector form to forecast the 
influence of multiple scattering effects on polarization characteristics of 
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the li&ar light signal. The method of statistical simulation (Monte-Carlo 
technique) forms the "basis of the algorithm. 

Numerical estimates obtained for a model of stratocumulus at 
X = 0.69^3]i under boundary conditions close to the conditions of natural 
experiment being discussed proved to be in a good agreement with the results 
of observation. 

Specifically, Fig. 1 shows the profiles of polarization (P) versus 
depth (L) of the following drop formations: fog (curve l) with horizontal 

meteorological visibility of U km two stratocumulus at a height of 1100 m 

—1 —1 

with the attenuation factors 6 = 0.01m (curve 2) and 6 = 0.05m (curve 3). 
Curve 3 shows the results of numerical estimates and the value of their 
statistical error. 



Fig. / 
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STATISTICAL PROPERTIES OF BISTATIC LIDAR SIGNALS 


J. A. Reagan, B. M. Herman, and D, J. MacKinnon 
The University of Arizona 
Tucson, Arizona 85721 
ABSTRACT 

Bistatic lidar systems which employ short pulse Q-switched lasers for 
transmitters yield "instantaneous” received signals that are quite analogous to 
those collected by monostatic lidar . In particular, the bistatic lidar 
equation for the instantaneous received power for scattering through an angle 
0 s is proportional to an instantaneous scattering depth, V , (depth along either 
the transmitter or receiver axes) given by 



where & is the laser pulse length. Bistatic signals are only received during 

the time the transmitted laser pulse is within the bistatic common volume, 

and' the common volume depth, L , along the transmitter axis is given by 

0 


2 r 2 6 B 
sin 0 

s 


where r^ is the distance from the center of the common volume to the receiver 
and 0^ is the half-angle of the conical receiver field-of-view. As the received 
bistatic signals will be independent from one another for scattering from in- 
stantaneous scattering volumes separated by approximately £' , the effective 
number of independent signal samples, N, collected during the passage of the 
laser pulse through the common volume is given by 

h : L t = 11 r a 6 b Blng (r-) . 

V % sin 0 

s 


If the rms signal to noise ratio of the instantaneous bi static signal is n. 
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time integration of the bistatic signal would therefore be expected to yield an 
improved signal to noise ratio of n /f . Assuming temporal stability, the 
signal to noise ratio may be further improved by averaging the integrated 

signals collected for a number of laser shots. 

\ 2 ] 

For the University of Arizona bistatic lidar system , the observed in- 
stantaneous rms signal to noise ratio is less than 5 for typical measurement 
conditions (9 from 100-150° and scattering heights 1-2 km), and theoretical 
calculations using representative aerosol distributions yield similar results. 
In order to improve the signal to noise ratio to a level which will permit 
meaningful aerosol size distribution information to be extracted from bistatic 
measurements , signal integration and averaging techniques of the type described 
above are required. Sample observations and supporting statistical calculation 
are presented to demonstrate the degree of accuracy one can practically attain 
with such signal improvement methods . 


1. Reagan, J. A., and B. M. Herman: Proc. ibth Radar Met. Conf. , 

pp. 275-280, (1970). ' - 

2. Reagan, J. A., and W. P. Webster: "A Bistatic Lidar for Measuring 

Atmospheric Aerosol Distributions Paper presented at 3rd Conference on 
Laser Radar Studies of the Atmosphere, Ocho Rios, Jamaica, W.I., (1970). 
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AEROSOL SIZE DISTRIBUTIONS DETERMINED FROM 
BISTATIC LIDAR OBSERVATIONS 

B. Herman 5 J. Reagan , T. King, R. Peck 
ABSTRACT 

Bistatic lidar observations at several scattering angles have been 
collected with the University of Arizona bistatic lidar (l). These 
'measurements are estimated to have errors not exceeding ± 3 % , thus 
representing a considerable improvement over earlier results reported 
with this system (2). These improved measurements have permitted 
mathematical inversions to be performed for determination of the aerosol 
size distribution (3). Results of such inversions will be presented. 

The technique for the mathematical inversions used in this work is 
departure from that employed in the earlier studies in that a modified 
form of the technique reported by Chahine (4) is used. This technique 
has the advantage of guaranteeing a nonnegative solution (i.e., the 
particulate number density is always positive within any size interval, 
as it must be) even in the presence of large observation errors. Of 
course the quality of the solution deteriorates as the errors increase, 
but nevertheless, an optimal solution results, A more complete description 
of this technique will be presented. 


Ill 
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BOUNDS ON THE IMAGINARY PART OF THE AVERAGE 
INDEX OF REFRACTION OF TROPOSPHERIC AEROSOLS 
FROM BISTATIC LASER SCATTERING 

By 

Gray Ward 

Department of Electrical Engineering 
and 

Department of Physics and Astronomy 
University of Florida 
Gainesville , Florida 32601 

ABSTRACT 

Measurements of the degree of polarization and the intensity as a function 
of scattering angle, compared with a theoretical model combining a regularized 
Junge distribution with Mie and Rayleigh scattering cross sections, has allowed 
us to determine the complex average index of refraction -of the atmospheric 
aerosol. The values of the real part of the index of refraction cluster around 
1.50, and lie within +.05. The values of the imaginary part of the index of 
refraction average -0.005, within a factor of two. Three parameters of the 
particle size distribution were also determined. Data taken during March 
through June 1972 will be presented along with best fits from an atlas of 
theoretical curves. Inaccuracies arise due to an oversimplified particle size 
distribution, and from the nonspherical particle problem. 
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MEASUREMENT OF RAINFALL INTENSITY BY LIDAR 


S. T. Shipley, E. W. Eloranta, J. A. Weinman 
Department of Meteorology, University of Wisconsin 
1225 West Dayton Street 
Madison, Wisconsin 53706 
ABSTRACT 

The University of Wisconsin monostatic lidar has been fired horizontally 
in rainfall, providing digital information on returned power in 15 meter in- 
crements to about 5 km. Simultaneous operation of rain gages under the beam 
path indicates that derived optical extinction and gage rainfall rates are 
correlated. Results from a stratus shower of 20 October, 1972 are given in 

figure I; the experimental relationship between the extinction coefficient 

-1 

3 (km ) and rainfall rate R (mm/ hr) is 

3 ~ 0.13 R' 9 ^ 

A relaxation algorithm giving the extinction coefficient from lidar returns 
as a function of range has been developed. Examples of lidar returns and the 
derived extinction coefficients from spatially inhomogenous showers are pre- 
sented. Experimental correlations of the extinction coefficient to rainfall 
rate are also given for spatially inhomogenous rainfall. 


Ilk 






Figure I. Average extinction as measured by lidar plotted 
with average rainfall rates given by three tipping bucket 
rain gages. The variance of the regression coefficients is 
due to the scatter of the data points about the line of 
regression. The correlation coefficient is 0.94. 
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MEASUREMENTS OF TRANSMISSIVITY OF LASER LIGHT THROUGH CLOUDS 

Ch. Werner 

German Establishment for Aeronautical and Space Research (DFVLR) 

Institute for Atmospheric Physics 
D-8031 Oberpfaffenhofen, Germany 
ABSTRACT 

By lidar measurements of clouds we determined that the transmissivity is 
better than it could be expected theoretically for the direct laser beam. For 
example, we-' 1 could measure a cirrus cloud through a thick altostratus layer. 

If the lower cloud was absent , the signal from the higher cloud was only two 
times greater. 

The scattering in clouds almost is a forward scattering. A portion of the 
scattered light has the same direction as the direct beam. This portion will 
be estimated for an empirical solution. The forward scattered portion of the 
extincted radiation in a layer d will be named a . The sum of directed 
(x) and diffuse transmissivity a(l - x) then is 

d 

x’ = x + a(l - x) , x = exp { - / adx } . 

0 

The Figure shows the result. The oscillograph photo of a lidar measurement is 

on the right side. By calculations with the lidar equation at first we get the 
2 

value B = 1 (fr = backscatter coefficient). For different portions 

2 

a(a = 0, , . . ,0, 5 > • • • j0.95 ) then we get the backscatter coefficient $ = B/t 1 
Figure 1 shows that an agreement with the measurement for the value a greater 
0.7 can be accepted. The standard visibility (V^ z ~ — ) is a 

value for the accuracy (upper scale). Visibilities of some hundred meters with- 
in the clouds may be correct . 

With other examples and model calculations the possibilities of utilization 
the forward scattering for range finding in fog and the limits for the other 
remote sensing methods are discussed. 
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AN ANALYSIS OF THE RISE TIMES 
OF LIDAR RETURNS FROM CLOUDS 


Ly 

John C. Alishouse 

National Oceanic and Atmospheric Administration 
National Environmental Satellite Service 
Suitland, Maryland 
and 

Sherman K. Poult ney 
Department of Physics and Astronomy 
University of Maryland 
College Park, Maryland 

ABSTRACT 

Analysis of lidar returns from clouds shows that the rise time of the 
lidar return varies as the cosecant of the lidar' s elevation angle. Plass 
and Kattawar have calculated the returned pulse shape for cloud models 
having both constant particle densities and monotonically increasing particle 
densities. Comparisons of the experimental values with the theoretical 
values shows the experimental values lying between the theoretical values. 
Calculations show that the increase in rise time attributable to the lidar 's 
beamwidth and changes in elevation angle are small in comparison with the 
observed changes in rise time. The observed changes in rise time could 
result from the presence of a boundary or developmental layer at the cloud 
base. 
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RAMAN-LAS'ER STUDIES OF THE LOWER ATMOSPHERE 

by. 

D. Ebeling, F. Frungel, G. Kuper 
Impulsphysik GmbH, Hamburg 

ABSTRACT 

The development of an air-pollution measuring system considering Raman- 
and Resonant Raman Scattering is under progress at different groups in 
Germany since 1971. The preliminary experiements are undertaken with a 
high-power high repetition rate ruby-laser system either at the uv or the 
red ruby line. The receiver consists of a conventional Cassegrain optical . 
system, followed either by vary narrow line filters (2-6 X) or double mono- 
chromators. Optimization of emitter and receiver optics, including photon- 
enhancement by TIR (total internal reflection light guides), showed same 
surprising results for the received Photon flux and SNR. An overall effi- 
ciency calculation is given for Ruby, N^ and Dye. Laser Systems, taking into 
account normal and Re sonant -Raman Scattering.' 

The commercially available long life time Ruby Laser was developed for 
other purposes at Impulsphysik GmbH, Hamburg, (150 MW, 20 ns, 2 cps), 
the Ng Laser at the University of Darmstadt and the Dye Laser at the ISL 
Laboratories of Weil/Rhein. Power ranges for both the Ng and Dye Laser 
are 0.1 - 10 MW at repetition rates of 10 - 100 cps. 

Using high gain low noise multipliers emission values in the range of 50 
100 ppm are to be measured by an A-scope technique and imission values of 0.2 
5 ppm with photon counting methods . The insertion of a multi-channel-plate 
imate intensifier - vidicon data processing system is taken into consideration 
as a future step. 

Quantitative measurements are carried out in an "artificial atmosphere" 
produced in defined ppm quantities in a large vacuum-vessel located some 
distance away from the experimental set up. 

Laser heads are constructed to be easily exchanged for different ware 
lengths. These measurements are running at Impulsphysik GmbH in Hamburg until 
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August 1973 in cooperation with the MBB laboratories and the University of 
Munich. In parallel exact scattering cross-section measurements are 
carried out at the Spectroscopic Institute of the University of Munich under 
the leadership of Prof. Dr. Brandmuller . 
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DETERMINATION OF TRANSPORT PARAMETERS IN THE ATMOSPHERIC 
BOUNDARY LAYER BY DOPPLER OPTICAL RADAR 

G. Benedetti-Michelangeli , CNR, Frascati, Italy 

F. Congeduti , University of Florence, Italy 

G. Fiocco, University of Florence, Italy 

ABSTRACT 

With the Doppler Optical Radar it is possible to measure, simultaneously, 
and at different heights in the lower troposphere, the vertical component of 
wind velocity, w , and the aerosol concentration n. 

Thus it is possible to separate the data into average and fluctuating 
components 


w = w + w' 
n = n + n' 

and therefore obtain the eddy flux 

(j> = n'w' 

Also, since the gradient d n is also measurable, it is possible to obtain 

8 z 

the diffusion coefficient of aerosols K 

The variable character of <j> can be studied as well as its spectrum. 
Results of recent experiment will be described. 
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